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COORDINATED  SCIENCE  LABORATORY 
ABSTRACT  OF 
PROGRESS  REPORT  FOR  SEPTEMBER,  OCTOBER,  NOVEMBER,  1962 

1.  Electric  Vacuum  Gyroscope 

In  the  course  of  improving  the  electric  vacuum  gyro,  it  was  found  that 
some  beryllium  rotors  outgassed  sufficient  nitrogen  to  render  them  unusable. 
Continued  gyro  performance  tests  are  being  conducted  with  deliberately  dis- 
turbed gyro  parameters  in  order  to  test  theoretical  performance  criteria. 
The  status  of  the  multi-electrode  gyro  system  is  described. 

H.  W.  Knoebel 

2.  Surface  and  Atomic  Physics 

A  technique  based  on  the  work  by  Mitchell  is  being  used  to  grow  thin 

57 
single  crystals  of  AgCl  doped  with  Co   .   An  rf-discharge  Cs  lamp  is  developed 

for  use  in  optical  pumping  experiments. 

E.  Luscher 
3  .    Computers 

CSX-1  operation  statistics  and  preparation  of  library  routines  are  reported. 
Connection  of  the  CSX-1  to  the  CDC  1604  equipments  is  completed.   The  CDC  1604 
operation  statistics  and  the  modifications  to  the  FORTRAN  and  CODAP  systems  are 
given.   Preliminary  studies  on  an  artificial  intelligence  program  are  reported. 

R .  M .  Brown 

4.  System  Theory 

Various  results  related  to  stability  and  control  optimization  were  obtained. 
An  algorithmic  synthesis  procedure  for  sequential  circuits  was  tried  on  the  CDC 
1604.   Several  configurations  of  self-repairing  systems  were  studied.   Necessary 
and  sufficient  conditions  for  the  realizability  of  terminal  capacity  matrices  for 
communication  nets  were  derived.   In  addition  to  eight  journal  articles,  four 
technical  reports  were  issued. 

J.  B.  Cruz,  Jr 

5 .  PLATO 

During  this  quarter  work  continued  in  testing  existing  educational 
material  in  the  areas  of  computer  programming  and  electrical  engineering. 
New  material  was  prepared  with  the  cooperation  of  the  Mathematics  Project, 


University  High  School,  in  the  area  of  modern  high  school  mathematics. 
Both  the  equipment  and  the  computer  program  are  nearly  completed  for  the 
PLATO  III  system.   This  system  is  designed  to  automatically  control  as  many 
as  thirty-two  student  sites  by  means  of  our  CDC  1604  computer.   Two  of  the 
PLATO  III  sites  are  expected  to  be  operational  early  next  quarter.   A  new 
computer  program  was  completed  which  allows  the  students  to  use  the  PLATO 
system  as  a  simulated  laboratory.   An  effort  was  started  this  quarter  to  record 
various  physiological  measurements  of  the  students  as  they  worked  through  the 
educational  material.   Some  data  have   been  recorded,  but  an  analysis  has  not 
been  completed. 

D.  L.  Bitzer 

6.  Vacuum  Instrumentation 

During  this  quarter  the  differential  ionization  gauge  has  been  sufficiently 
developed  to  demonstrate  its  potential  usefulness.   Evaluation  of  the  suppressor 

gauge  (Pat.  Pend.)  has  continued  this  quarter.   Anomalous  behavior  of  the 
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Bayard-Alpert  gauge  under  certain  conditions  at  pressures  in  the  10    Torr 

range  is  being  investigated.   Comparative  studies  of  total  and  partial  pres- 
sures using  two  types  of  mass  spectrometers  have  continued. 

W.  C.  Prothe 

7.  Plasma  Physics 

In  the  course  of  assembling  instrumentation  for  the  experiment  on  elastic 
collisions  of  electrons  with  excited  atoms,  a  Polarity  Coincidence  Detector 

has  been  constructed  and  tested.   It  works  in  the  range  20  c/sec  -  20  Kc/sec 

-4 
and  is  capable  of  detecting  signal-to-noise  ratios  as  low  as  10   .   The  cross 

section  for  elastic  collision  of  electrons  for  Xenon  (Xe)  has  been  determined 

o 
in  the  temperature  range  of  200-600  K.   The  1604  program  for  the  solution  of 

the  Boltzmann  equation  for  shock  waves  has  been  completed  and  new  data  for  a 

Machnumber  of  3.17  has  been  obtained.   For  the  plasma  delay  line,  conditions 

for  the  constancy  of  the  time  of  flight  in  cylindrical  geometry  are  discussed. 

M.  Raether 

8.  High  Magnetic  Field  Superconductors 

In  further  studies  of  the  superconducting  solenoid  it  has  become  apparent 
that  there  are  serious  drawbacks  to  the  design  chosen.  A  program  is  in  pro- 
gress to  prepare  crystalline  Nb  Sn .   Work  is  in  progress  to  study  the  effect 


of  structural  defects  on  the  superconducting  properties  of  Nb  and  V.   A 
liquid  helium  level  indicator  has  been  perfected  and  submitted  for  patent. 
A  cooperative  program  has  been  instituted  with  Fansteel  Metallurgical  Corp. 
for  preparation  and  evaluation  of  some  high  field  superconductors  . 

C.  Satterthwaite 

9.    Vacuum  Electrical  Breakdown 


Electron  microscopy  with  improved  resolution  allows  the  observation  of 
points  on  electrode  surfaces,  both  for  contact  welds  that  have  been  separated, 
and  for  smoothing  of  cathode  points  upon  heating.   For  electrodes  in  ultra 
high  vacuum,  pre-breakdown  pulses  and  pre-breakdown  electrode  erosion  is 
observed.   Current-voltage  relationships  are  being  studied. 


E  .  Lyman 


Electron  emission  through  thin  (300  Angstroms)  films  of  Al  0  ,  under 


R.  N.  Peacock 


10.   Physical  Electronics 
Electron  emission  th 
high  stress  (10-15  volts),  is  reported 
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Systems,  IRE  Trans,  on  Electronic  Computers,  vol.  EC-11,  Aug.  1962,  pp.  459- 
465. 
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M.  Sobral,  Jr.,  "Stability  Considerations  in  Synthesis  of  Time-Varying 
Feedback  Systems,"  accepted  for  publication  in  the  AIEE  Transactions. 

J.  Resh,  "The  Inverse  of  a  Nonsingular  Submatrix  of  an  Incidence  Matrix," 
accepted  for  publication  in  the  IRE  TRANS,  on  Circuit  Theory. 

2.    Meeting  Papers  Presented 

D.  Bitzer  and  P.  Braunfeld,  "Description  and  Use  of  a  Computer  Controlled 
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1 1 lino is  . 
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1.1   Introduction 


ELECTRIC  VACUUM  GYROSCOPE 


H.  W.  Knoebel  R.  W.  Anderson  V.  Krone 

D.  0.  Skaperdas  J.  D.  Gooch  D.  Lee 

D.  Allen  B.  D.  Kirkwood  W.  C.  Prothe 


The  electric  vacuum  gyroscope  is  a  high  precision  inertial  navigation 
instrument  employing  a  spherical  metal  rotor  spinning  without  physical 
contact  while  supported  in  ultrahigh  vacuum  by  servo-controlled  electric 
fields.   The  rotor,  initially  accelerated  by  rotating  magnetic  fields, 
coasts  during  extended  periods  of  gyro  operation.   The  gimballed  housing 
technique  is  employed  to  minimize  spurious  torques  which  would  otherwise 
result  from  imperfect  rotor  sphericity.   Optical  data  readout  from  a  zig- 
zag pattern  on  the  rotor  is  employed  to  maintain  the  gimballed  housing  in 
accurate  alignment  with  the  rotor  spin  axis. 

Since  demonstration  of  the  electric  vacuum  gyroscope  feasibility  has  been 
accomplished,  the  new  objectives  of  this  group  are:   (1)  to  obtain  a  complete 
understanding  of  various  subtle  effects  pertinent  to  the  design  of  ultrahigh 
precision  electric  vacuum  gyroscopes  and  (2)  to  develop  new  techniques  for 
exploitation  of  the  electric  vacuum  gyroscope  principle. 
1.2   New  Gyro  Assembly 

In  the  previous  progress  report  it  was  mentioned  that  a  second  version 
of  the  gyro,  in  which  calibrated  displacements  and  misalignments  can  easily 
be  introduced,  was  being  readied  for  use.   Three  beryllium  rotors  obtained 
from  American  Beryllium  Company  of  Sarasota,  Florida,  were  to  be  used  to  test 
this  gyro.   Two  of  these  were  balanced  by  Speedring  Corporation  of  Warren, 
Michigan.   The  third  is  being  kept  in  reserve.   These  rotors  were  fabricated 
by  joining  two  halves  in  a  helium  atmosphere,  using  a  small  amopnt  of  MMM  C-4 

epoxy  in  the  joint.   All  three  rotors  passed  CSL's  helium  leak  test;  i.e., 

—8 
no  helium  leak  could  be  detected  using  a  sensitivity  of  10   Torr  liter/sec. 

Sinusoidal  readout  patterns  were  sand  blasted  on  the  two  balanced  rotors,  and 


one  of  the  rotors  was  placed  in  a  gyro  housing  for  pumping  down.   It  was  found 

— fi 
that  the  gyro  housing  pressure  could  not  be  reduced  below  about  10   Torr  with 

a  25  liter/sec  Varian  Vaclon  pump  of  2.5  liter/sec  effective  pumping  speed. 


When  the  Vaclon  pump  was  removed  from  the  housing,  leaving  only  the  combina- 
tion titanium  ion  pump  and  Bayard-Alpert  gauge  connection,  the  pressure  in- 

-4 
creased  to  about  10   Torr  even  after  numerous  flashings. 

After  re-examining  both  the  gyro  rotor  and  housing  to  determine  that  no 
external  leaks  had  developed,  the  first  rotor  was  replaced  by  the  second  one 
with  identical  results.   The  housing,  electrodes,  and  ceramic  spacers  were  also 
replaced,  with  identical  results,  leaving  us  to  suspect  that  the  rotors  were 
the  source  of  gaseous  contamination.   Accordingly,  one  of  the  rotors  was  sub- 
jected to  a  mass  spectrometer  analysis  using  a  Davis-Vanderslice  spectrometer. 

Spectrographs  were  taken  under  three  conditions:   (1)  of  the  rotor  before  a 

o  o 

bakeout,  (2)  of  the  rotor  after  a  bakeout  of  100  C.  for  18  hours  and  150  C. 

for  38  hours,  and  (3)  of  the  spectrometer  system  after  the  rotor  was  removed. 

Under  the  first  two  conditions  large  amounts  of  mass  28  (N  ,  CO  or  both),  N, 

and  A  were  observed,  the  relative  amount  of  mass  28  being  0.98,  that  of  N 

-3  -4 

being  3.5  x  10  ,  and  that  of  A  being  9.6  x  10  ,  with  a  total  pressure  of 

5.6  x  10   Torr.   Under  condition  (3)  no  N  or  A  was  observed,  while  the  rela- 
tive amount  of  mass  28  (only  CO  in  this  case)  was  about  .65,  with  a  total  pressure 

-8 

of  2.5  x  10   Torr.   Hence,  we  can  conclude  that  large  amounts  of  N  and  N   and 

some  A  are  degassing  from  the  rotor  in  such  a  manner  that  the  bakeout  described 
above  has  no  influence.   The  second  rotor  appears  to  behave  similarly.   Both 
rotors  arced  badly  when  we  tried  levitating  them.   The  third  rotor  held  in 
reserve  will  be  tested  for  vacuum  compatability  before  it  is  balanced.   Studies 
are  under  way  to  establish  the  reasons  for  the  behavior  of  this  particular 
batch  of  rotors  in  an  effort  to  avoid  future  repetition.   Previous  rotor 
fabrications  and  balancings  by  the  same  companies  showed  no  such  anomalous 
behavior . 

D.  Skaperdas 
1.3   Orientation  Sensing  Data  Processor 

A  packaged  version  of  the  superheterodyne  Orientation  Sensing  Data  pro- 
cessor was  put  into  operation  during  the  last  quarter.   For  improved  perfor- 
mance and  circuit  simplification,  some  changes  were  incorporated  into  this 
system.   The  conventional  pre-amplif iers  were  replaced  by  automatic  gain 


A  description  oi  the  first  system  is  given  in  C.S.L.  Progress  Report  for 
December  1961,  January  and  February  1962,  p.  17. 


compressing  pre-araplif iers  to  compensate  for  variations  in  the  photomicro- 
scope  output  .   The  eight-diode  phase  detector  in  the  phase-locked  reference 
loop  was  replaced  with  an  inverted-transistor  phase  detector.   The  frequency 
converter  in  the  reference  loop  has  been  changed  from  the  type  7360  beam 
deflection  tube  circuit  to  the  less  complicated  inverted-transistor  circuit. 

An  effort  was  also  made  to  simplify  the  converter  circuitry  in  the  in- 
formation channels.   However,,  in  this  case,  the  desired  information  signal 
is  approximately  60  db  smaller  in  amplitude  than  the  large  second-harmonic 
component  of  the  pattern  waveform.   Therefore,  the  double-balancfed  circuitry, 
as  used  in  the  7360  converters,  is  desirable  in  order  to  obtain  additional 
rejection  of  the  large  second-harmonic  component,  to  improve  the  signal-to- 
noise  ratio,  and  to  avoid  overloading  the  filters  that  follow  the  converters. 
For  these  reasons,  the  converter  circuitry  was  not  modified  in  the  information 
channels  . 

H.  V.  Krone 
1 .4  Gyro  Test  Data 

Continued  tests  are  being  conducted  to  determine  the  effects  on  the 
locus  of  the  gyro  spin  axis,  of  constant  east  or  west  translation  of  the 
centrifugally  distorted  rotor,  with  respect  to  the  housing.   Theory  predicts 
that,  for  this  type  of  rotor  perturbation,  the  spin  axis  locus  should  show 
spiraling  and  culmination  shifts,  with  the  closure  time  not  affected.   In 
all  of  these  runs  the  data  agreed  with  theory.   The  theoretical  results  are 
shown  in  Figure  1  in  which  spiraling  rate  and  culmination  for  a  given  value 
of  u   are  shown  versus  the  nominal  angle  at  which  the  spin  axis  cones  about 
north.   The  coefficient  w   represents  the  maximum  drift  rate,  in  radians  per 
day,  of  the  rotor  spin  axis  due  to  a  constant  displacement  S   of  the  center  of 
mass  from  the  center  of  electric  support,  where 


w    =   S.mg/iw  . 
si       1     r 

The  value  of  S   depends  on  the  amount  of  centrifugal  deformation  of  the  rotor 

and  the  gap  space  between  the  rotor  and  electrodes .   The  curves  of  Figure  1 
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Figure  1.   Gyro  spiraling  and  culmination  shift 


o 
for  coning  angles  greater  than  90  .   The  coefficient  co   has  been  calculated 

and  tabulated  for  hollow,  thin  shells  and  solid  aluminum  and  beryllium  rotors 

for  various  rotor  speeds  in  CSL  Report  M-95 .   The  curves  for  beryllium  rotors 

are  reporduced  in  Figure  2,  where  Ax/d   in  that  figure  corresponds  to  SA/d»  and 

Ap  corresponds  to  u   in  minutes  of  arc  per  day . 

To  illustrate  the  use  of  these  curves,  consider  an  actual  run  in  which  a 

hollow  beryllium  rotor,  with  a  preferred  moment  of  inertia  and  spinning  at 

240  rps,  was  used.  This  rotor  was  given  a  constant  eastward  displacement, 

S  ,  with  respect  to  its  electrodes,  of  700  microinches  out  of  a  gap  spacing 

of  7  mils  .   Referring  to  Figure  2,  and  choosing,  as  a  fair  approximation, 

that  curve  lying  midway  between  a  hollow  and  solid  sphere  (since  our  rotor 

is  thicker  about  its  equator),  then,  for  Ax/d  =  0.1,  we  obtain  oo   ~  70 

arc  min/day  =  0.02  rad/day.   If  the  gyro  spin  axis  is  given  an  initial  con- 

o 
dition  such  that  its  locus  comes  about  north  with  a  nominal  radius  of  15  , 

then,  from  Figure  1,  the  spiraling  rate  for  oo  -   70  arc  min  per  day  will  be 
340  x  0.02  =6.8  arc  min  per  day,  and  the  culmination  shift  will  be  420  x  .02 
=8.4  arc  min.   The  spiral  direction  and  culmination  shift  will  be  in  the  direc- 
tions indicated  in  the  upper  left  hand  corner  of  Figure  1.   The  actual  data 
agreed  very  closely  with  the  above  theoretical  results. 

The  curves  of  Figure  1  were  obtained  from  the  theoretical  Equations  (10) 
and  (11)  in  CSL  Progress  Report  for  December  1961,  January  and  February  1962 
for  our  actual  latitude  L  =   40.1   instead  of  45  ,  and  with  the  following  cor- 
rection.  In  Equation  (7)  of  that  report,  the  torque  due  to  rotor  displacements 
S   and  S   should  be 

Tb  =  R  X  (mg.Q)  Q, 


while  that  due  to  displacement  S  ,  as  before,  is  identical  to  the  mass  un- 
balance torque.   With  this  change,  and  for  a  general  latitude  L,  the  spin 
axis  locus  equations  become 


— —  =  -  (P  Sin  L  -  B  Cos  L)  [go  +  oj  „  -  w  (P  Cos  L  +  B  Sin  L)  ] 
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si  s2 


in  which 


A      =        Ml  -   a2) 


B      = 


vli  -  q2  -  p2) 


D .    Skaperdas 


1.5  Digital  Computer  Data  Analysis 

It  was  mentioned  in  the  preceding  Progress  Report  that  the  coefficients 
of  a  theoretical  set  of  data  could  not  be  exactly  reproduced  by  the  computer, 
when  the  set  of  data  was  fed  to  the  Computer.   The  program  used  is  described 
in  the  preceding  Progress  Report  .   The  difference  between  the  actual  data  and 
the  theoretical  curve  generated  by  the  computer  contains  a  periodic  variation. 
Hence,  either  the  program  is  not  achieving  a  proper  fit  or  Equations  (1.1)  and 
(1.2)  are  not  accounting  for  all  torques  acting  on  the  gyro. 

Approximate  solutions  in  closed  form  of  Equations  (1.1)  and  (1.2),  using 
first  order  perturbation  theory,  have  been  obtained.   For  conditions  under 
which  our  gyro  is  normally  operated,  these  solutions  agree  very  closely  with 
the  numerically  obtained  exact  solutions  of  Equations  (1-1)  and  (1.2).   A 
computer  program  is  currently  being  written  to  use  these  new  solutions  to 
obtain  the  six  coefficients  in  Equations  (1.1)  and  (1.2)  from  actual  data. 

D.  Lee,  D.  Skaperdas 

1.6  Multi-Electrode  Strapped-Down  Gyro  Analysis 

The  obvious  advantages  of  a  gyro  housing  being  mounted  firmly  to  its 
carrier  has  led  to  an  investigation  of  the  possible  advantages  of  a  multi- 
electrode  housing  for  this  application.   Relatively  large  precession  torques 
may  appear  if  the  spin  axis  of  a  distorted  rotor  is  free  to  move  relative  to 
the  electrodes.   A  gimballed  housing,  slaved  to  the  spin  axis,  minimizes 
these  torques  . 


Most  undesirable  precession  torques  occur  because  the  rotor  is  not  a 
perfect  sphere,  and  electric  field  forces  do  not  all  converge  at  the  rotor 
center.   The  shape  of  the  rotor  with  cylindrical  symmetry  about  the  spin  axis 
may  be  expressed  as  the  Fourier  series 

oo 

a  =  a   +  2  a   cos  n6. 
o     „  n       ; 
n=l 

where  a  is  the  rotor  radius,  a   is  the  undistorted  radius,  and  0  is  the  co- 

'      o  ' 

latitude  angle  measured  from  the  spin  axis  .   The  a   are  coefficients  des- 

n 

cribing  the  magnitude  of  each  of  the  n  distortions,  which  are  divided  into 
two  classes,  even-order  and  odd-order  distortions  (n  even  and  odd,  respec- 
tively) .   An  analysis  is  being  made  of  the  precession  torques  which  exist  on 
a  distorted  rotor  placed  in  a  non-gimballed  multi-electrode  housing.   Some  of 
the  conclusions  which  have  been  drawn  from  this  analysis  follow. 

In  a  non-gimballed  electrode  housing,  a  fundamental  requirement  for 
minimum  precession  torque  is  that  the  first-order  force  distribution  obtains 
over  the  electrode  sphere.   The  first-order  force  distributions  are  defined 
in  the  Progress  Report  for  December  1961,  January  and  February  1962.   A 
necessary  condition  for  this  distribution  is  that  the  sum  of  forces  exerted 
on  the  rotor  by  each  electrode  pair  remains  constant,  independent  of  external 
forces,  and  is  the  same  constant  for  all  electrode  pairs.   Electrodes  must 
be  arranged  over  the  sphere  in  opposing  pairs  and  each  electrode  of  a  pair 
must  be  identical  in  shape.   If  this  condition  is  met,  plus  two  other  condi- 
tions, that  the  rotor  be  perfectly  centered  in  the  electrode  sphere,  and 
that  the  effect  of  the  electrode  boundaries  be  negligible,  then  precession 
torques  on  a  rotor  with  only  even-order  distortions  are  nulled.   Under  these 
same  conditions,  precession  torques  on  a  rotor  with  only  odd-order  distor- 
tion do  not  null,  but  consist  of  (1)  first-order  precession  torques  which 
add  to  the  mass  unbalance  torque,  and  (2)  higher-order  torques  which  decrease 
approximately  as  the  inverse  square  of  the  total  number  of  electrodes  . 

Finite  interelectrode  spacings,  which  are  necessary  to  prevent  voltage 
breakdown,  introduce  additional  torques  on  a  rotor  with  both  even  and  odd 
order  distortion.   These  torques  are  approximately  proportional  to  the  pro- 
duct of  the  effective  spacing  and  the  electrode  diameter. 


The  precession  torques  which  occur  on  a  non-centered  rotor  with  even- 
order  distortion  are  approximately  the  same  as  torques  occurring  on  a  per- 
fectly centered  rotor  with  odd-order  distortions .   These  torques  decrease  as 
the  inverse  square  of  the  electrode  number. 

To  verify  the  foregoing  conclusions,  an  experiment  is  being  set  up  in- 
volving a  purposely  distorted  rotor  placed  in  both  a  housing  with  six  round 
electrodes,  and  then  a  housing  with  twenty  triangular  electrodes,  the  latter 
placed  with  minimum  electrode  spacing  according  to  the  pattern  of  an  icosa- 
hedron .   The  electric  torques  will  be  determined  by  noting  the  equilibrium 
position  of  the  "spin  axis"  of  the  purposely  distorted  non-spinning  rotor  as 
a  function  of  the  orientation  of  that  axis  with  respect  to  the  movable  elec- 
trode housing.   The  electric  torques  are  counterbalanced  by  the  known  mass- 
unbalance  torque  which  is  the  rotor's  first-order  distortion. 

The  six-electrode  housing  is  the  same  assembly  which  is  presently  used 
in  CSL's  gyro  program.   Figure  3  is  a  photograph  of  the  twenty-electrode 
housing,  along  with  the  toroidal   inductors,  which  are  an  integral  part  of 
the  resonant  support  servo  mechanisms  to  be  used  in  this  experiment. 

D.  Allen 
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2.1  Growth  of  AgCl-Single  Crystal 

2 
A  technique  based  on  that  of  Clark  and  Mitchell   is  being  used  to  grow 

thin  discs  of  single  crystals  of  AgCl,  doped  to  \1°   with  Fe  CI  .   These  single 

crystals  are  to  be  used  in  a  Mossbauer  experiment,  which  will  study  the  iron 

in  the  AgCl  environment . 

The  crystals  are  grown  from  commercial  AgCl  powder  in  a  "furnace  consisting 
of  two  stainless  steel  slabs,  of  dimensions  3  in.  x  3  in.  x  1/4  in.  and  12  in.  x 
3  in .  x  1/4  in.,  placed  edge-to-edge  separated  by  a  1/16  in.  gap.   The  smaller 
slab  and  the  adjacent  end  of  the  larger  slab  are  placed  on  separate  heating 
coils,  the  entire  apparatus  being  placed  inside  a  pyrex  cylindrical  glass 
enclosure  4  in.  o.d.,  which  sits  on  a  transite  base.   The  ends  of  the  glass 
enclosure  are  sealed  off  with  two  "clamping"  endplates,  made  of  stainless 
steel,  with  0- rings  and  nuts  and  bolts  doing  the  clamping.   Two  needle  valves 
exit  from  these  plates,  permitting  the  growing  of  the  crystals  in  a  CI   gas 
atmosphere . 

The  crystal  itself  is  grown  on  a  quartz  glass  disc,  of  1.5  in.  diameter. 
After  the  powder  is  placed  on  this  disc  and  melted,  a  second  thin  pre-heated 
quartz  disc  is  placed  on  top  of  the  AgCl,  being  separated  from  the  first  disc 
by  three  thin  glass  rods .   This  results  in  a  thin  AgCl  melt  being  formed 
between  the  two  quartz  discs  . 

The  discs  are  then  placed  on  the  smaller  steel  slab  which  is  preheated 

o  o 

to  500  C.   The  near  end  of  the  larger  slab  is  heated  to  350  C,  resulting 

in  a  negative  temperature  gradient  down  the  larger  slab,  the  far  end  being 

o 
at  100  C.   The  discs  are  drawn  slowly  toward  and  down  the  length  of  the 

larger  slab,  using  an  external  motor  and  pulley  system,  and  a  connecting 

steel  rod   arm  .   As  the  melt  crosses  the  gap,  a  crystal  forms  and  grows 


2 
Journal  Photographic  Science,  4,  1-20,  (1960) 
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o 
on  the  cooler  side  of  the  gap,  the  melting  point  of  AgCl  being  450  C.   The 

discs,  upon  reaching  the  100  C  end  after  the  5  hour  drawing  time,  are  quickly 

immersed  in  distilled  water,  freeing  the  crystal  from  the  glass.   At  present, 

some  AgCl  single  crystals  have  been  obtained. 


2.2   High  Intensity  Cs-Lamp  for  Optical  Pumping 

A  high  intensity,  low  noise  cesium  light  source  has  been  developed  for  use 

3 
in  optical  pumping  studies.   The  lamp  design  proposed  by  Brewer   for  rubidium 

has  been  modified,  and  many  experimental  bulbs  have  been  prepared  and  tested 

in  order  to  produce  maximum  intensity  in  the  D   resonance  line  of  cesium. 

The  total  radiant  power  has  been  measured  to  be  greater  than  4.2  watts  under 

normal  operating  conditions  .   This  may  be  compared  with  a  widely  used  cesium 

4 
light  source  produced  commercially  by  the  Varian  Company  which,  at  maximum 

5 

power,  develops  a  total  radiant  power  of  0.14  watts  . 

The  power  in  the  single  D   cesium  line  is  of  particular  interest  in 

optical  pumping.   We  have  incorporated  the  present  lamp  into  a  previously 

6 
described  apparatus  ,  and  have  measured  the  optical  pumping  signal  for  the 

cesium  Zeeman  resonances.   We  compare  this  signal  to  those  obtained  using  a 

7 
Varian  lamp,  and  an  Osram   lamp.   The  geometry  was  optimized  in  each  case, 

and  for  the  Varian  lamp  a  special  reflector  was  used.   The  results  are  sum- 
marized in  Table  I .   Details  concerning  the  construction  and  operation  of  this 
new  cesium  lamp  are  included  in  a  paper  now  in  preparation,  which  will  be 
issued  as  a  CSL  report . 


3 

R.  G.  Brewer;  Rev.  Sci .  Instr.  32,  1356,  (1961) 

4 
Research  Cesium  Spectral  Lamp  available  from  Varian  Associates,  611  Hansen  Way, 

Palo  Alto,  California. 

5 

All  measurements  of  light  flux  were  made  using  a  calibrated  Eppley  thermopile. 

The  intensity  from  the  Varian  lamp  was  made  at  the  maximum  power  suggested  by 
the  manufacturer. 

6  •' 

W.  W.  Holloway,  E.  Luscher,  R.  Novick;  Phys  .  Rev.  126,  2109,  (1962). 

7 
Available  from  Edmund  Scientific  Co.,  Barrington,  N.  J. 
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Table  I .   Comparison  of  Various  Cesium  Lamps 


Source 


Optical  Pumping  Signal 
(Arbitrary  Units) 


Total  Radiant 
Power  (watts) 


Source  Area 
(cm  ) 


Varian 
Osram 
Present  Work 


1.5 

15 
100 


0.14 


4.2 


3 
20 
25 


14 
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COMPUTER 


L.  Kypta  S.  Srivastava 

J.  Leichner  J.  Stifle 

M.  Levitt  R.  Trogdon 

W.  W.  Lichtenberger  F.  Wise 

P.  Ponzo  N.  Wax 


3.1   Introduction 

This  division  is  concerned  with  the  design  and  application  of  digital 
computers  for  information  processing,  primarily  in  the  areas  of  real  time  and 
non-numerical  operations.   The  present  work  lies  in  three  areas:   the  develop- 
ment and  supervision  of  the  CSL  computing  facilities  centering  on  the  CSX-1 
and  a  CDC  1604  computer,  and  studies  in  artificial  intelligence. 

3  .2   CSX-1  Computer 

3.2.1  Operations 

Period:   September  1  to  November  29,  1962 

Total  Running  Time:  714  hrs . 

Ave.  per  day  (5  day  week):  10.5  hrs. 

Scheduled  Maintenance  Time:  0.0  hrs. 

Emergency  Maintenance  Time:  8.5  hrs.  (1.2%) 

Operational  Time:  705.5       (98.8%) 

Breakdown  of  Emergency  Maintenance  Time 

Paper  Tape  Reader:  7  hrs. 

Power  Supply  1.5  hrs. 

3.2.2  Input-Output 

The  1612  line  printer  and  the  1607  magnetic  tape  units  have  been 
connected  to  the  CSX-1.   A  Flexowriter  has  also  been  connected  to  the  CSX-1, 
providing  an  additional  input-output  facility.   All  of  these  units  have  been 
checked  out  and  are  fully  operational. 

J.  Stifle 
R.  Trogdon 
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3.2.3   Programming 

Programming  efforts  for  the  CSX-1  during  the  last  quarter  were  directed 
toward  planning  the  future  SAC  II  assembly  system.   The  assembling  principles  are 
to   be   changed  in  order  to  utilize  the  equipment  to  be  made  available  to 
the  computer.   Furthermore,  portions  of  the  assembly  program  will  be  rewritten 
to  accomodate  future  programming  requests  and  to  take  advantage  of  the  proposed 
extensions  in  the  order  code. 

Now  completed  and  available  in  SAC  I  are  library  routines  to  read  and 
punch  data  in  Flexowriter  code,  and  to  perform  useful  binary-to-Flex  code 
translation.   Library  subroutines  to  read  and  write  magnetic  tape,  to  output 
to  the  1612  line  printer,  and  to  perform  other  necessary  code  translation 
are  now  being  developed. 

Two  programs  for  the  CSX-1  were  also  finished  during  this  period.   The 
first  is  a  general  integer  -  fraction  printing  routine  described  in  the  CSX-1 

Programming  Supplement  I  .   The  second  program  is  a  random  number  generator 

8  th 

using  the  following  algorithm  for  the  (i  +  1)    random  number: 

N   ,  =  (2?  +  1)  N.  +  5,168,620,107  x  2~31 

l+l  i         '         ' 

N   =  0  or  can  be  specified  by  the  user.   A  31  bit  number  (two  CSX-1 

registers)  is  generated.   Successful  tests  were  run  on  the  most  significant 

16  bits  of  the  numbers  generated.   The  tests  used  were  a  sexadecimal  frequency 

distribution  test,  serial  correlation  test  with  lags  of  1  to  15,  runs  of  l's 

and  0's  vertically  through  successive  numbers  at  each  digit  position,  and  runs 

horizontally  through  successive  words  . 

R .  M .  Brown 
R.  D.  Jenks 
F.  Wise 

3.3   CDC  1604  Computer 

3.3.1   Operations 

Period:   September  1  to  November  30,  1962 

Total  Running  Time:  1036    hrs . 

Average  per  day  (five  day  week)  16.44  hrs. 

Average  per  day  (including  Saturday,  Sunday, 

and  holidays)  11.38  hrs. 


A.  Rotenberg,  "A  New  Pseudo-Random  Number  Generator,"  JACM  7,  pp.  57-77. 
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Operational  Time:  955.17  hrs . 

Preventive  Maintenance  Time:  79    hrs. 

Emergency  Maintenance  Time:  1.83  hrs. 

Actual  Operating  Time/Scheduled  Operating  Time:  99.79  % 

Breakdown  of  Emergency  Maintenance  Time: 

Main  frame  of  computer:  1.5   hrs. 

Punch  power  supply:  0.33  hrs. 

E.  H.  Neff 

3.3.2   1604  Programming 

In  the  past  quarter  some  minor  difficulties  in  the  FORTRAN  system 
have  been  detected.   These  difficulties,  enumerated  below,  have  all  been 
corrected,  and  the  new  systems  tape  will  be  installed  on  January  2,  1963. 
The  trouble  dump  routine  was  observed  to  give  incorrect  results  when  dumping 
high  areas  of  core.   Routine  bindump  was  found  to  place  extraneous  informa- 
tion in  core  under  certain  circumstances.   These  errors  were  corrected. 

It  has  been  the  practice  to  use  the  systems  tapes  with  the  file  protect 
rings  removed.   It  was  observed  that  routine  TRANSFER  made  use  of  the  write 
backspace  and  so  could  not  be  used  conveniently.   The  routine  was  changed 
to  make  use  of  the  read  backspace . 

The  operation  of  that  portion  of  the  resident  routine  which  reads  cards 
(or  their  images)  was  found  not  to  be  consistent  with  its  description  in  the 
manual.   Appropriate  corrections  were  made. 

Finally,  a  need  has  been  felt  for  some  time  for  a  user  to  be  able  to 
print  out  relative  time  when  operating  under  monitor  control.   A  routine 
for  doing  this  has  been  written. 

An  executive  system  which  will  enable  the  CSX-1  to  be  used  as  an  on- 
line satellite  of  the  1604  has  been  worked  out  and  many  of  the  routines  coded, 
Final  assembly  and  de-bugging  remain.   This  project  is  being  done  mostly  as 
an  experiment  in  techniques  of  time-sharing  and  computer  intercommunication. 
It  is  not  currently  planned  to  use  the  on-line  satellite  other  than  for  a 
short  period  of  evaluation. 

The  Illinois  Systems  Tape  has  been  in  use  for  one  quarter  now  without 
incident  .   Some  minor  revisions  and  improvements  have  been  made  to  two  of 
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the  routines  in  ILLRES,  the  utility  package  which  accompanies  the  system. 
The  SNAPSHOT  dump  routine  has  been  extended  so  that  it  now  has  the  ability 
to  print  out  fixed-point  fractions,  integers,  and  normalized  floating- 
point numbers  in  addition  to  its  previous  capabilities.   Work  is  nearly  com- 
pleted on  a  change  in  program  LIST  which  will  produce  a  dictionary  of  sym- 
bolics and  their  references . 

W.  Lichtenberger 
3.4  Artificial  Intelligence 

During  this  quarter  a  series  of  programs  was  written  to  investigate 
computer  learning  techniques.   One  of  the  problems  in  this  area  is  the  assign- 
ment of  credit  (or  negative  credit)  when  a  long  series  of  decisions  results  in 
a  correct  (or  incorrect)  solution.   The  approach  here  is  to  assume  that  the 
influence  of  decisions  near  the  end  of  the  decision  tree  can  be  determined 
with  greater  confidence  than  decisions  further  removed  from  the  final  state. 
When  these  'terminal"  decisions  have  been  evaluated  with  high  confidence 
(as  the  result  of  successive  approaches  to  the  problem)  then  the  decisions 
which  led  to  the  terminal  decisions  can  be  evaluated  by  the  same  process. 

The  particular  game  used  as  a  vehicle  for  this  study  is  an  as  yet 
unnamed  variation  of  Nim .   Two  players  successively  remove  markers  from  an 
N  x  N  array,  a  legal  move  being  the  removal  of  one  or  more  markers  from  a 
row  or  column.   The  player  removing  the  last  marker  loses  the  game. 

There  were  several  reasons  for  choosing  this  game.   One  reason  was  that 
while  the  rules  of  play  are  simple,  it  is  certainly  not  a  trivial  game  if 
N  >  4.   No  algorithm  of  correct  play  is  known  to  exist.   Moreover,  the  game 
is  basically  one  of  geometric  patterns  and  under  the  rules  of  play  many 
patterns  can  be  shown  to  be  equivalent  under  rotation,  translation  and  inter- 
change of  rows  (or  columns).   This  means  that  a  wealth  of  "primitive'  patterns 
(65,536  for  a  4  x  4  array)  can  be  reduced  to  a  much  smaller  number  of  canonical 
forms  (201  for  the  same  array).   This,  of  course,  simplifies  the  treatment  of 
the  game  in  the  computer. 

In  preliminary  experiments  of  pitting  the  computer  against  human  opponents 
of  varying  skill,  the  machine  did  demonstrate  learning.   Generally,  about  thirty 
games  were  sufficient  to  bring  machine  performance  to  the  point  where  it  would 
win  most  games . 

J.  Divilbiss 


4.   SYSTEM  THEORY 
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4.1  Introduction 

This  group  is  concerned  with  three  general  areas  of  research.   One  area 
is  automatic  control  theory.   For  the  most  part,  attention  is  focused  on 
various  aspects  of  optimum  and  adaptive  control  systems  .   The  synthesis  of 
stable  control  systems  where  various  types  of  performance  indexes  are  opti- 
mized is  a  major  objective.   Another  area  of  research  is  circuit  theory  in- 
cluding linear  graph  theory  and  computer  methods  in  network  design.   Applications 
of  linear  graph  theory  to  communication  nets  and  switching  theory  are  being 
studied.   A  third  area  of  research  is  sequential  machine  theory.   Presently, 
work  is  in  progress  on  algorithmic  design  of  sequential  switching  circuits, 
survey  of  automata  theory,  and  self -repair. 

J.  B.  Cruz,  Jr 

4.2  A  Method  of  Control  Law  Generation  Using  Sequential  Phase  Space  Constraints 
During  the  quarter,  R-154  and  R-155  were  written  and  published.   This 

marks  the  end  of  the  first  phase  of  this  work. 

Several  approaches  to  the  computational  aspects  of  the  optimization 
problem  have  been  considered  and  a  few  modifications  of  existing  methods  have 
been  examined  in  some  detail.   At  the  present  time  no  conclusive  results  can 
be  reported. 

Work  during  the  next  quarter  will  entail  detailed  investigation  of  a 
performance  criterion  modification  technique  recently  proposed  and  an  attempt 
to  implement  it  computationally.   Application  of  the  adaptive  constrained  descent 
technique  to  optimization  problems  will  be  continued  and  numerical  examples 
will  be  solved. 

S.J.  Kahne 
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4.3   Use  of  Lyapunov  Functions  in  the  Synthesis  of  Stable  Optimum 

Linear  Systems 

If  the  performance  of  a  linear  discrete  system  can  be  adequately  evaluated 
by  a  quadratic  function  of  the  state  variables,  and  if  the  constraints  on  the 
allowable  control  vector  can  be  expressed  as  saturation-type  constraints,  then 
a  synthesis  procedure  can  be  given  that  results  in  stable  optimum  design  for  a 
very  wide  class  of  sampled  data  control  problems.   The  synthesis  technique 
described  here  results  in  a  feedback  control  that  is  relatively  easy  to  calcu- 
late and  simple  to  implement,  in  marked  contrast  to  results  obtained  using 
dynamic  programming  or  the  maximum  principle.   A  technical  report  presenting 
these  results  in  detail  is  in  preparation.   Only  the  highlights  of  the  solu- 
tion will  be  given  here. 

We  first  formulate  the  problem  quantitatively.   The  discrete  plant  is 
represented  by  the  state  equations 

x(n+l)  =  cb(x)  x(n)  +  A(X>  u(n),  (4.1) 

~        Tn   ~       n       ' 

(x) 

where  x  is  the  state  vector,  u  the  control  vector,  <f>    the  transition  matrix, 

(x)  ~  n 

and  A    the  distribution  matrix.   Moreover,  the  system  outputs  are  related  to 
n 

the  states  by  the  equation 

c(n)  =  M(X)  x(n).  (4.2) 


We  assume  the  system  input  process  may  be  represented  as  the  output  of  a  dis- 
crete dynamic  system: 

w(n+l)  =  4,(W)  w(n)  +  A(W)  s(n).  (4.3) 

~        ^n   ~       n   ~  * 

(w) 

where  w  is  a  state  vector,  s  a  process  input  vector,  (f>    a  transition  matrix, 

(w)  ~  n 

and  A    a  distribution  matrix.   The  output  of  this  system,  r(n),  is  related 

n  ~ 


to  the  state  by 


(w) 
r(n)  =  IVT  '   w(n),  (4.4) 


20 


We  can  combine  these  equations  into  a  single  set  of  equations  by  defining 


4>  = 

Jn 


(x)   „ 

d>      0 

•  n 


M  = 


M(X)    0 


0      M 


(w) 


J 


A, 


(x)' 


0    J 


A 


(w) 


(4.5) 


z(n)  = 


x(n) 


w(n) 


2<n> 


c(n) 
r(n) 


(4.6) 


Then  the  equations  become 


z(n+l)  =  H>  z(n)  +  A  u(n)  +  ^  s(n), 

~        Yn  ~       n  ~       n  ~ 


(4.7) 


y(n)  =  M  z(n) 


The  control  vector  u(n)  is  subject  to  constraints: 


I  u  (n)  I  <  c., 
J      -J 


(4.8) 


for  all  n. 

The  performance  measure  is 


Y=  S   [c(n)-r(n)]  Q '  [  r(n)  -  c(n)  ], 
n=l  n 


(4.9) 


where  Q'  >  0.   We  may  write  this  in  terms  of  the  state  vector  z(n)  as 
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J  =  2  z  (n)  Q   z(n), 
n=l 


(4.10) 


where 


(M(X))T 


(-M(W))T 


Q'  [M(X)   -M(W) 


(4.11) 


We  may  now  state  the  problem  concisely  as  follows:   For  the  system 


z(n+l)  =  *  z(n)  +  A  u(n)  +  ^  s(n), 

~        Tn  ~       n  ~       n  ~    ' 


(4.12) 


y(n)  =  M  z(n), 


choose   u(n)    subject    to   constraints 


u .(n)   I     <   c. 
J  -      J 


(4.13) 


so  as  to  minimize 


J  =  2  z  (n)  Q   z(n), 

n  ~ 
n=l 


(4.14) 


and  so  that  the  resulting  feedback  system  is  stable. 

The  author  has  obtained  the  solution  to  this  problem  for  a  variety  of 
cases.   The  system  may  be  time-varying.   The  inputs  may  be  deterministic  or 
random.   The  system  parameters  may  be  varying  in  a  random  fashion.   Finally, 
the  state  may  not  be  known  exactly,  but  must  be  estimated  from  noisy  data. 
For  the  random  cases,  the  expected  value  of  Equation  (4.9)  is  minimized,  and 
the  resulting  system  is  stable  in  the  mean.   The  method  makes  use  of  Lyapunov 
functions  to  yield  a  system  that  is  at  once  stable  and  optimum.   The  details 
of  the  technique  are  too  lengthy  to  present  here,  and  will  be  the  subject  of 
a  forthcoming  technical  report.   In  all  cases  mentioned  the  control  turns  out 
to  be  of  the  form 
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u  .  =  -c  .  sat  [C  z  ]  .,  (4.15) 

J      J        ~  J 


where 

+  1,     for   [C  z]  j>  1, 

sat  [C  z].  :=  I    [C  z]  ,    fori  [C  z]  j  I  <  +  1,  (4.16) 

for    [C  z].   <  -1 
~  J 

The  state  variable  feedbacks  to  the  sat-function  are  thus  merely  simple  linear 
gains  .   If  z  must  be  estimated,  it  turns  out  that  the  estimation  may  be  carried 
out  separately.   The  control  for  the  overall  optimum  system  still  being  given 
by  Equation  (4.12)  with  z  replaced  by  the  estimated  z. 

Future  investigation  will  involve  possible  extensions  of  this  technique 
to  the  control  of  nonlinear  systems.   In  connection  with  this  contemplated 
extension,  work  has  been  done  on  methods  of  generating  Lyapunov  functions  for 
nonlinear  systems.   Certain  methods  related  to  energy  (stored  energy,  total 
energy,  Hamiltonians,  etc.)  seem  promising. 

W.  R.  Perkins 
4.4   Closed  Loop  Frequency  Response  of  Systems  with  More  than  One  Nonlinear 

Element 

A  number  of  methods  have  been  proposed  so  far  for  obtaining  the  closed 

loop  frequency  response  of  nonlinear  systems.   Almost  all  the  methods  proposed 

9, 10, 11 

use  the  concept  of  describing  functions  to  represent  the  nonlinear  elements.  '   ' 

Also,  invariably,  nonlinear  systems  with  only  one  nonlinear  element  has  been 
considered.  (See  Figure  4.) 


9„ 

A  Frequency  Response  Method  for  Analyzing  and  Synthesizing  Contractor  Servo- 
mechanisms".  R.J.  Kochenburger,  AIEE  Transactions,  Vol.  69,  Part  I,  (1950). 

"A  Generalized  Method  for  Determining  the  Closed  Loop  Frequency  Response  of 
Nonlinear  Systems,"  L.T.  Prince,  Jr.,  AIEE  Transactions,  Part  II,  Vol.  73, 
September,  1954. 

Closed  Loop  Response  of  Nonlinear  Systems  by  a  Functional  Transformation 
Approach,"  J  .C .  Hill,  IRE  Transactions,  PGAC,  Vol.  AC-7,  Number  4,  July  1962. 
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The  closed  loop  "transfer  function"  for  a  given  sinusoidal  input  R(co) 
is,  in  general, 

_,  ,      N(E,co)  GT (u) 

C(co)  L  , 

wri  = '  (4-17) 

(<JJ)         1  +  N(E,w)  G  (co) 

Li 

where  the  describing  function  of  the  nonlinear  element  is  a  function  of 
both  the  amplitude  and  frequency  of  its  input, 

|£4  =  N(E,oo).  (4.18) 

E(u>)    - 

The  open  loop  transfer  function  of  the  system  is 

C(oo) 


E(oo)       '  L 

From  Equations  (4.17)  and  (4.19)  we  get 

E(u>) 


=  N(E,co)  G  (oo)  .  (4.19) 


R(oj)    1  +  N(E,w)  GT  (co) 

Li 


(4.20) 


To  obtain  the  closed  loop  frequency  response,  evaluation  of  the  ratio 


C(u) 


at  constant  amplitude  R  and  different  frequencies  is  required,  whereas  only 

C(w) 

—-. — r-  at  constant  error  amplitude  is  known,   Since  N  is  a  function  of  E  and  go, 

C  (co) 
a  direct  calculation  of  the  ratio  — — -  is  not  possible.   A  graphical  construe- 

R(w)  12 

tion,  which  is  one  of  the  methods  proposed  so  far,    can  be  used  however. 

This,  in  summary,  involves  plotting  families  of  curves.   One  is 


C 

E 


E  constant 


for  different  values  of  E  from  Equation  (4.19),   Curves  obtained  are  superim- 
posed on  Ni chol's  chart  and 


12 
""Analysis  of  Nonlinear  Control  Systems,"  D.  Graham  and  D.  McRuer,  John  Wiley 

and  Sons,  1961 . 
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C 
R 


E  constant 


curves  are  obtained.   By  cross  plotting  the  families  of  these  curves  I  —  I 
versus  go  for  R  constant  can  be  obtained. 

The  application  or  extension  of  this  method  to  systems  with  more  than  one 
nonlinear  element  fails,  because  the  closed  loop  transfer  function  of  such  a 
system  involves  products  of  describing  functions.   For  example,  in  Figure  5, 
one  has 


N   =  N  (E,u), 

(4.21) 
N2  =  N2(C,gg), 


rf    S  N  (E,oj)  G  (u) 

C(oo)  1       L  (4.22) 


R(go)    1  +  N  (E,go)  N  (C,oj)  G  (go) 
1        2        1j 


The  feasibility  of  using  an  iteration  procedure  is  considered  in  this  note. 
The  investigation  has  not  been  completed.  The  outline  of  this  procedure  is 
as  follows: 

For  the  sake  of  simplicity  it  is  assumed  that  N   and  N   are  frequency 
invariant,  i.e.,  N   =  N  (e)  and  N   =  N  (C) . 
Step  1 

In  Figure  5,  assume  that  the  feedback  path  is  not  closed.   For  a  given 

value  of  R(oj)  obtain  the  output  C(go)  =  C  (go)  .   This  can  be  done  because  N   = 

o  1 

N  (R  )  is  known . 
Step  2 


Compute  N  using  C   as  input.   Assume  feedback  is  now  closed.   Compute 

r(    ^        Ni(R)  Gt  (co> 
C(go)  1     L 


R(go)    1  +  N(R)  N(C  )  GT  (go) 
1     2   o   L 


For  the  given  value  of  R,  compute 
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3, 
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or 
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C(w)  =  C1(oj) 


Step  3 

Using  C  (u)  as  input  to  N  ,  compute  N   =  N  (C  ) ,  and 
X  2t  £  £•       J. 


r<    ^  VR)  Gt(cj) 

C(oj)  1      L 


R(u>)  "  1  +  N1(R)  N2(C1)  GL(W)  * 


As  long  as  the  feedback  signal  is  small,  it  is  assumed  that  N  remains  un- 
affected.  Continue  the  iteration  procedure  until  successive  values  of  C(co) 
are  reasonably  close  to  each  other.   Repeat  for  different  frequencies  to 
obtain  the  frequency  response  characteristic, 

As  mentioned  earlier,  this  investigation  has  not  been  completed.   While 
preliminary  verification  by  numerical  examples,  involving  small  feedback 
signals,  point  to  the  convergence  of  the  iteration  procedure,  the  validity 
of  this  method  with  large  feedback  signals  remains  to  be  checked. 

V.  Gourishankar 
4.5   Systems  Optimization  as  a  Game 

The  use  of  game  theory  to  solve  systems  optimization  problems  have  been 
suggested  by  several  authors.   '   '     Kalman    regards  the  optimization  of 
a  linear  second  order  stationary  system  as  two-person,  zero-sum  game  between 
Nature  (which  attempts  to  maximize  the  minimum  transient  time  by  choosing 
the  initial  state  of  the  system)  and  Man  (who  tries  to  minimize  the  maximum 
transient  time  by  choosing  the  damping  factor)  .   He  shows  that  the  opponents 
have  good  pure  strategies  and  opens  the  question  whether  this  is  true  in  games 
of  this  type . 


13 

R.  E.  Kalman,  "Control  System  Analysis  and  Design  Via  the  'Secona  Method'  of 

Lyapunov,"  Journal  of  Basic  Engineering,  A.S.M.E.,  June,  1960,  pp.  371-393. 

14  1  11  "B 

J  .G .  Truxal  and  J.  J.  Padalino,   Decision  Theory,   Chapter  15  of   Adaptive 

Control  Systems,"  edited  by  E.  Mishkin  and  L,  Braun,  Jr.,  McGraww-Hill  Book 

Co.,  New  York,  (1961)  . 

15 

R.  Bellman,  "Adaptive  Control  Processes:   A  Guided  Tour,"  Princeton  Univ , 

Press,  Princeton,  New  Jersey,  (1961)  . 

16 

Reference  13,  cited  above,  pp .  388-389. 
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17 
Bellman    also  considers  Nature  as  an  intelligent  opponent  and  using 

the  principle  of  optimality  finds  a  computational  algorithm  to  solve  the  problem. 

14 
Truxal  and  Padalino    point  out  the  application  of  decision  theory  in  the 

evaluation  and  design  of  feedback  systems .   In  doing  so,  they  consider  feedback 

problems  and  communication  problems  as  a  statistical  game  (that  is  a  game  in 

which  one  of  the  opponents,  Nature,  is  not  considered  as  a  conscious  opponent 

18 
who  can  take  advantage  of  mistakes  made  by  the  other  opponent   )  .   Some  recent 

19 
research  in  this  field  has  been  done  by  Horing.     In  this  report  the  Bang- 
Bang  Control  problem  (restricted  to  linear  plants)  is  treated. 

Consider  now  the  following  optimization  problem.   Given  the  system  state 
equation 

x  =  f  (x,u,v),  (4.23) 

where  x  is  a  nxl  state  vector,  u  is  the  control  vector  (mxD  and  v  is  the 
disturbance  vector  (pxl)  • 

Let  u  be  restricted  to  a  certain  vector  space  U,  i.e.,  u  is  constrained 
to  be  a  vector  in  U.  Let  v,  be  restricted  to  a  certain  vector  space  V,  i.e., 
v  e    V.   Let  the  performance  criterion  be 

t 

M(t1)  =       g[x(t),u(t),v(t)  ]  dt.  (4.24) 

t 
o 

Assume  that  the  purpose  of  the  design  is  to  minimize  Equation  (4.24)  by  choosing 
an  appropriate  u(t) . 


17 

R.  Bellman,  "Adaptive  Control  Processes:   A  Guided  Tour,"  Princeton  Univ. 

Press,  Princeton,  New  Jersey,  (1961),  pp.  180-193. 

18 

D.  Blackwell  and  M.A.  Girshick,  "Theory  of  Games  and  Statistical  Decisions," 

John  Wiley  and  Sons,  Inc.,  New  York,  (1954). 

19 

S.  Horing,  "On  the  Optimization  Design  of  Predictor  Control  Systems,"  Research 

Report  No.  PIB  MRI-945-61,  Polytechnic  Institute  of  Brooklyn,  Microwave  Research 

Institute,  August,  1961. 
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The  problem  just  formulated  can  be  considered  as  a  statistical  game, 
as  follows.   The  game  G  is  a  triple  G  =  (V,U,M)  where  player  I  (Nature) 
chooses  v  £  V,  player  II  (the  engineer)  chooses  u  £  U  and  the  pay-off  is 
M(t  )  . 

In  the  next  quarter,  this  approach  to  systems  optimization  is  to  be  under- 
taken.  First,  simple  cases  for  instances,  restricting  the  systems  to  be  linear, 
time  invariant,  and  U  to  be  the  set  of  vectors  u(t)  such  that  u(t)  -  1,  or 
u(t)  =  -1,  or  u(t)  =  0,  will  be  considered  in  order  to  achieve  a  better  under- 
standing of  the  problem,   Also  Kalman's  question  will  be  carefully  studied, 

M,  Sobral,  Jr, 
4.6   Design  of  Linear,  Stable  Controllers  for  Linear  Multivariabie  Systems 

4.6.1  Introduction 

The  problem  of  designing  linear,  stable,  physically  realizable  controllers 
for  continuous  or  discrete  linear  stable  multivariabie  plants  subject  to  deter- 
ministic or  stochastic  inputs,  in  order  to  satisfy  a  certain  performance  index 
with  a  certain  constraint  is  being  investigated „ 

The  methods  available  at  present  for  the  synthesis  of  the  controllers  do 
not  give  stable  controllers  in  general.   Therefore,  we  may  have  to  use  analog 
computers  with  positive  feedback  to  realize  unstable  controllers  obtained  by 
the  existing  methods.   It  is  proposed  to  design  stable  controllers  which  in 
general  can  be  realized  at  lower  costs.   Purther,  the  constraint  imposed  in 
this  investigation  has  not  been  considered  in  previous  works. 

The  configuration  in  Figure  6  is  assumed  for  the  system  where  P  is  the  n- 
input,  n-output  fixed  plant  and  G  and  H  are  the  controllers  to  be  designed. 
The  following  discussion  is  concerned  with  continuous  systems  with  stochastic 
inputs.   The  methods  can  be  extended  for  discrete  systems  with  deterministic 
inputs  . 

4.6.2  Statement  of  the  Problem 

Given  the  plant  transfer  function  matrix  P,  the  input  correlation  matrix 

$   ,  the  input-desired-output  cross-correlation  matrix  $    and  the  cesired 
rr'  rd' 

output  correlation  matrix  $  ,,  we  have  to  design  linear,  stable,  physically 

T 
realizable  controllers  G  and  H  such  that  Av[e  (t)L  e(t)]  is  minimized  where 

6(t)  =  C(t)  -  D(t),  D(t)  is  the  desired  output,  and  C   the  output  due  to 

~       ~R       ~      ~  ~N 

noise  is  kept  below  a  certain  level.   Here  P,  $   ,  $  , ,  and  $,  are  n  x  n 

'   rr'   rd'       dd 
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matrices  and  are  functions  of  the  variable  S  where  S  =  jw.   Here  we  have  assumed 

the  existence  of  the  Fourier  Transforms  of  P,  $   ,  $   and  3?,,.   L   is  a  n  x  n 

'   rr    rd      dd     1 

positive  definite  matrix  of  constants  and  C(t)  and  D(t)  are  n  x  1  matrices.   We 
have  assumed  that  the  noise  is  acting  at  the  output  terminals  of  the  plant. 
The  noise  entering  anywhere  in  the  system  can  be  replaced  by  an  equivalent 
noise  N  at  the  output  of  the  plant  as  we  have  done  here. 
Referring  to  Figure  6,  we  write 


=  (I+PGH)~   N  +  (I+PGH)"1  PG  R, 
=  W  N  +  WPG  R, 
where  W  =  (I+PGH)"  . 


We  notice  here  that  the  ideal  solution  for  W  would  be  a  zero  matrix  in 
order  to  get  complete  rejection  of  noise  in  the  output  .   This  will  lead  to 
unrealizable  controllers.   Let  W  be  any  specified  matrix  different  from  the 
zero  rejection  matrix,  in  order  to  get  a  required  rejection  of  noise  in  the 
output  . 

4.6.3   Restrictions  on  the  Specification  of  W 

The  system  in  Figure  6  can  be  reduced  to  the  system  in  Figure  7  as  far 
as  the  signal  transmission  is  concerned. 

Since  the  system  is  to  be  stable,  P  is  stable,  and  G  is  designed  to  be 
stable,  either  W  should  be  stable,  or  the  controller  G  should  have  as  its 
zeros  the  unstable  poles  of  W .   A  matrix  is  said  to  be  stable  if  all  its 
elements  are  stable.   By  saying  that  G  should  have  as  its  zeros  the  unstable 
poles  of  W,  we  mean  that  every  element  of  matrix  G  should  contain  as  its  zeros 
the  unstable  poles  of  the  elements  of  W .   If  W  is  unstable,  the  output  due 
to  noise  will  increase  with  time,  and  hence  we  assume  W  to  be  stable. 

In 

(I  +  PGH)  =  W_1, 

we  want  G  and  H  to  be  stable  controllers,  and  since  P  is  stable,  W    should 
be  stable.   Solving  the  above  equation  for  H,  we  have 
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H  =  G-1  p  1    (W_1  -I). 

If  H  is  to  be  stable  W    -  I  should  contain  as  its  zeros  the  R*H-P-  zeros  of 
det  P  and  det  G. 
Summarizing, 

(1)  W  should  be  stable, 

(2)  W~   should  be  stable, 

(3)  W  '  -  I  should  contain  as  its  zeros  the  R-H-P-  zeros  of  det  P,  and 

(4)  W  "  -  I  should  contain  as  its  zeros  the  R-H°P°  zeros  of  det  G. 
W  can  be  specified  to  satisfy  conditions  1,  2,  and  3.   However,  this  W 

may  not  satisfy  condition  4.   It  seems  that  W  cannot  be  specified  to  satisfy 
condition  4,  since  G  is  not  known,  and  to  design  G,  W  has  to  be  specified  to 
start  with.   If  condition  4  is  not  satisfied,  H  is  not  stable.   If  the  con- 
troller H  is  not  necessarily  to  be  stable,  W  has  to  satisfy  only  condition  1. 

4.6 .4   Solution  for  the  Special  Case 

If  we  assume  W  and  P  commute,  i.e.,  WP  =  PW,  then  all  the  four  conditions 
can  be  easily  satisfied, 

In  Figure  8  a  cascaded  systems  is  shown  where  P  is  the  plant  and  G'  is  the 

T 

controller  to  be  designed  to  minimize  Av[e  (t)Le(t)]  where  e(t)  =  C   (t)  -  D(t). 

20  1  ~      ~R 

This  G'  can  be  solved  by  using  the  equation 

[$   G,T  PT  P]+  =  [«  J  P]+, 
rr  rd     ' 

where  P  =  P(-S)  and  [  ]   indicates  the  operation  of  making  a  partial  fraction 

expansion  in  terms  of  left  half  plane  poles  only.   The  G'  so  obtained  is  stable, 

since  we  have  assumed  the  existence  of  Fourier  Transforms  of  $    <$    and  P. 

rr'   rd 

The  overall  system  transfer  function  of  the  system  in  Figure  8  is  PG '  and  this 

system  gives  the  minimum  mean   square  error  for  the  specified  plant  P,  $  , 

and  <£     The  plant  in  Figure  8  is  modified  to  be  WP  -  PW  where  W  is  a  stable 
rd 

matrix.   The  controller  G  for  this  plant  WP  =  PW  as  in  Figure  7  and  for  the 

same  3>   and  $  ,  as  given,  will  be  W  G',  and  this  is  stable  since  W    is 
rr      rd  '  21 

specified  to  be  stable.   This  fact  can  be  seen  from  the  solution   of  the 
equations 


20 

R.C.  Amara,   The  Linear  Least  Squares  Synthesis  of  Continuous  and  Sampled 

Data  Multivariable  Systems,  '  Technical  Report  No.  40,  Stanford  Electronics 

Labs.,  P.  59;  Eq .  3-29. 

21 

Ibid,  p.  61,  Eq.  3-39. 
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[*   GT  (PW)  (PW)T]+  =  [$rd  PW]+ 


[*   G,TPT  P]+  =  [*  J  P]+, 
rr  rd 


The  overall  system  transfer  function  of  the  system  in  Figure  7 ,  which  is  the 
system  transfer  function  of  Figure  6  as  far  as  the  signal  transmission  is 
concerned,  is  PW  W   G'  =  P  G' .   This  transfer  function  is  independent  of  W, 
and  hence  the  mean  square  error  is  independent  of  W. 

Substituting  G  =  W  G'  in  the  solution  for  H  we  have 

H  =  G~1p~1(W_1  -  I) 

=  G,~1WP-1(W~1  -  I)  =  G'~1P~1(I-W);, 

Since  W  and  P    commute.   Since  G'  and  P  are  known,  even  before  W  is  specified, 
(I-W)  can  be  specified  to  contain  the  R-H-P-  zeros  of  det  P  and  det  G1. 

4.6 .5  Future  Work 

A  method  of  specifying  W  to  satisfy  condition  4  or  to  give  a  minimum 
phase  G  (that  is  det  G  does  not  have  any  R-H'P*  zeros)  is  being  studied.   Any 
specified  W  has  to  satisfy  all  the  four  conditions  in  order  to  design  stable 
controllers,  and  hence  a  given  W  should  be  modified  to  satisfy  all  the  four 
conditions.   Also,  the  effect  of  the  specification  of  W  on  the  mean  square 
error  is  being  studied. 

4.6.6  Design  of  Forward  Controller  with  Feedback  Controller  Specified 
In  the  above  procedure  it  was  assumed  that  (I  +  PGH)    =  W  is  a  speci- 
fied matrix.   However,  we  are  not  sure  that  the  above  matrix  is  the  best 
specification  for  W.   Therefore,  we  assume  some  statistical  information  about 

the  noise  and  find  a  stable,  physically  realizable  W  which  minimizes 

T  T 

Av[e  (t)L  €(t)]  +  Av[N  (t)L  C  (t)].   We  also  assume  that  the  controller  H 

is  specified  and  H  and  H    are  stable.   Here  L   and  L   are  positive  defi- 
nite matrices  of  constants,  and  e(t)  =  C  (t)  -  D(t),  where  D(t)  is  the 
desired  output  . 
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The  system  in  Figure  6  can  be  reduced  to  the  system  in  Figure  9,  where 

A  -    I-W.   Using  the  calculus  of  variation  techniques,  we  have  the  following 

T  T 

equation  for  A,  which  minimizes  Av[e  (t)L  e(t)]  +  Av[C  (t )L  C  (t ) ] : 


[H_1(-S)   $rr(S)  [H_1(S)  f   AT(S)  L^+ 

=  [H-1(-S)   $   (S)   L   +$   (S)   L]+. 
rd       1    nn      2 

Again  we  are  assuming  the  existence  of  the  Fourier  transforms  of  H,  H   , 

$   ,  and  $   ,  and,  hence,  the  solution  of  the  above  equation  for  A  should 

rr'       nn      '       ' 

be  stable  and  physically  realizable.   The  solution  for  A  will  be  in  terms  of 

L  ,  the  LaGrange  multiplier  matrix,  and  the  value  of  L   can  be  solved  by 

T 

letting  Av[C   (t)L  C  (t)]  <  a  be  some  specified  value. 

From  this  A,  we  have 

G  =  P_1  A  (I  -  A)"1  H~  . 

The  G  so  obtained  is  not  necessarily  stable,  because  P   and  (I-A)    are 
not  necessarily  stable. 

4  .6  .7   Future  Work 

The  effect  of  change  in  specification  of  H  on  A  is  being  studied,  and 

this  can  be  used  to  design  G  as  a  stable  controller. 

T 
A  new  formulation  for  solving  for  G  directly  to  minimize  Av[e  (t)L  e(t)]  + 

T 
Av[C„T  (t)L  C„(t)]  is  being  studied. 
N      2  N 

R.  Narayanasamy 
4.7   On  Oriented  Communication  Nets 

In  this  past  quarter  a  note  was  submitted  to  the  PGCT  which  will  appear 
in  the  March  1963  issue  of  the  Transactions .   The  note  establishes  a  close 
connection  between  the  inverse  of  the  incidence  matrix  of  an  oriented  tree 
and  a  certain  path  matrix  of  the  tree. 

The  problem  of  the  realizability  of  matrices  as  terminal  capacity  mat- 
rices of  oriented  communication  networks  is  still  under  investigation.   Nec- 
cessary  and  sufficient  conditions  for  the  realizability  of  a  large  class  of 
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asymmetric  matrices  have  been  obtained.   It  is  suspected  that  the  conditions 
are  in  fact  general,  but  this  has  not  been  established  yet . 

J.  A.  Resh 

4.8  Survey  of  Papers  on  Graph  Theory 

With  the  increasing  application  of  graph  theory  to  the  electrical  net- 
work, a  collection  of  papers  on  pure  graph  theory  is  important  .   Many  workers 
have  made  contributions  to  the  scientific  aspects  of  graph  theory.   Unfor- 
tunately, these  published  papers  are  scattered,  and  some  are  not  readily  ob- 
tainable.  However,  more  than  two  hundred  papers  concerning  this  subject  have 
been  collected,   One  of  the  major  objectives  in  collecting  these  papers  was 
to  make  them  readily  available  for  the  Circuit  Group.   These  papers  will  be 
read  subsequently,  and  the  possibility  of  further  application  to  electrical 
networks  will  be  investigated. 

Wai-Kai  Chen 

4.9  Realizability  Conditions  of  Single-Contact  (SC)  Functions 

If  a  given  function  F,  which  is  expressed  as  a  union  of  products  of 
variables,  is  an  SC  function,  it  must  satisfy  certain  conditions  which  are 
known  to  be  the  necessary  conditions  of  an  SC  function.   For  example,  it  is 
known  that  any  SC  function  must  satisfy  the  "odd-ring-sum"  condition:   the 
ring  sum  of  an  odd  number  of  product  terms  in  the  given  function  F  must  be 
either  a  term  of  F  or  a  redundant  term  which  contains  all  the  variables  in 
some  term  of  F .   It  is  also  known  that  a  function  which  satisfies  the  odd- 
ring-sum  condition  may  not  be  an  SC  function.   However,  if  a  given  function 
does  not  satisfy  any  of  the  known  realizability  conditions,  it  cannot  be  an 
SC  function.   Thus,  a  method  of  testing  for  the  non-realizability  of  a  given 
function  to  be  an  SC  function  can  be  developed  from  a  set  of  necessary  con- 
ditions for  SC  functions. 

Let  F  be  a  function  of  e  variables,  expressed  as  a  union  of  products  of 
the  e  variables  .   Then  the  following  statements  are  each  the  necessary  con- 
ditions of  an  SC  function: 

If  F  is  an  SC  Function  it  is  necessary  that 

(1)  F  satisfy  the  odd-ring-sum  condition, 

(2)  F  contain  at  least  two  disjoint  basic  cut-sets  if  it  does  not 
contain  any  terms  of  degree  1, 
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(3)   the  number  of  disjoint  basic  cut-sets  in  F  must  not  exceed  the 
least  degree  among  all  the  terms  of  F. 

If  F  satisfies  all  the  necessary  conditions  listed  above,  certain 
"functional  transformations"  may  be  applied  to  it  in  order  to   transform 
it  into  a  modified  form.   Through  such  transformation,  either  the  given  func- 
tion is  considerably  simplified  so  that  conventional  synthesis  methods  may 
be  applied  to  it,  or  it  is  transformed  into  such  a  form  that  the  non-realiz- 
ability  of  the  original  function  may  be  detected  with  the  help  of  some  simple 
tests.   Thus,  by  applying  the  tests  developed  from  the  above  necessary  con- 
ditions to  the  given  function,  before  resorting  to  the  conventional  synthesis 
methods,  considerable  amount  of  time  and  effort  may  be  saved. 

The  details  of  various  tests  in  the  testing  method,  as  well  as  the  func- 
tional transformations  mentioned  above,  may  be  found  in  the  report,  Realiza- 
bility  Conditions  of  Single-Contact  Functions,"  soon  to  be  available. 

S.  P.  Chan 

4.10  Linear  Graph  Theory 

There  are  many  practical  problems  related  to  planar  graphs.   However, 
the  only  property  of  such  graphs  known  today  is  the  existence  of  dual  graphs . 
The  search  for  further  properties  of  planar  graphs  is  underway. 

W  .  Mayeda 

4.11  Design  of  Sequential  Switching  Circuits 

Conventional  methods  for  the  design  of  sequential  switching  circuits 
become  too  cumbersome  even  for  moderately  large  circuits.   Work  is  in  pro- 
gress on  a  method  for  algorithmically  designing  circuits  from  the  input- 
output  specifications  .   A  machine  language  program  is  being  written  for  the 
CDC  1604  computer  to  verify  the  capabilities  of  the  method. 

J.  Elsey 

4.12  Automata  Theory 

a)    A  survey  of  work  done  on  the  field  was  undertaken  six  months  ago 
to  determine  potential  avenues  of  research  which  may  be  pursued  here.   A 
report  on  this  survey  together  with  an  annotated  bibliography  will  be  issued 
in  the  near  future . 

H.  Y.  Chang 
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b)    It  has  been  found  that  automata  which  contain  only  one  output 
state  have  some  interesting  properties.   It  has  been  possible  to  charac- 
terize the  sets  of  tapes  accepted  by  such  automata.   Further,  any  regular 
set  of  tapes  can  be  characterized  in  terms  of  the  special  class .   Further 
work  is  being  done  along  these  lines  with  a  view  to  producing  a  canonical 
form  for  a  regular  expression. 

S  .  Seshu 
4.13  Self-Repair 

The  task  of  evaluating  different  possible  organizations  of  a  self-re- 
pairing system  was  undertaken  in  the  present  quarter.   It  is  proposed  to  simu- 
late these  organizations  on  the  CDC  1604  computer,  and  evaluate  them  by  in- 
jecting random  failures. 

J  .  Kruus 
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R.  Blomme  J.  Easley  B.  W.  Voth 
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5.1  Introduction 

The  purpose  of  the  PLATO  project  is  to  develop  an  automatic  teaching 
system  capable  of  tutoring  simultaneously  a  large  number  of  students  in  a 
variety  of  subjects.   The  central  control  element  of  the  teaching  system 
is  a  general  purpose  digital  computer.   The  PLATO  system  differs  from  most 
teaching  systems  in  that  a  single  high  speed  digital  computer  is  used  to 
control  all  student  stations.   Thus,  it  has  available  the  power  of  a  large 
digital  computer  in  teaching  each  student. 

5.2  Introduction  to  Programming  ILL I AC 

During  this  quarter  19  students, who  were  enrolled  in  a  senior  mathe- 
matics course  at  University  High  School,  participated  in  the  second  phase 
of  the  study  using  the  seven  chapter  sequence,  'introduction  to  Computers." 

Eight  students  have  already  completed  the  full  seven  chapter  sequence 
and  the  data  from  their  lessons  are  being  prepared  for  analysis.   The  other 
eleven  students  are  at  various  stages  in  the  sequence  and  should  complete 
it  by  the  end  of  December.   About  50  hours  of  ILLIAC  time  have  been  used 
to  date  for  this  phase  of  the  experiment , 

Booklet  reproductions  of  the  chapters  are  being  distributed  upon  the 
completion  of  each  lesson,  and  the  students  have  often  remarked  that  they 
do  use  the  booklets  for  further  study  and  review.   If  the  student  has  been 
unable  to  finish  a  chapter  in  the  allotted  PLATO  time,  he  has  completed  it 
from  the  booklet  before  his  next  PLATO  lesson.   This  has  resulted  in  much 
less  use  of  ILLIAC  time  for  "review"  purposes. 

0.  R.  Brown,  Jr , 

5.3  Preparation  of  Materials  for  Use  in  Teaching  Training 

The  four  members  of  the  University  of  Illinois  Committee  on  School 
Mathematics  (UICSM)  staff  who  were  engaged  25%  time  with  CSL  this  quarter 
have  spent  their  time  learning  about  the  PLATO  system  and  selecting  topics 
for  putting  into  an  instructional  sequence. 


39 


Topic  choices  from  UICSM  materials  have  been  made  (e.g.  classifying 
expressions  via  given  pattern-forms)  and  the  instructional  sequence  is  in 
outline  form. 

0.  R.  Brown,  Jr. 

5.4  Lessons  in  Maxwell's  Equations 

The  three  lessons  on  an  Introduction  to  Maxwell's  Equations,  which  were 
prepared  last  year  by  Prof.  E.  C.  Jordan,  were  used  again  this  fall  for  seven 
students  from  Prof.  W.  G.  Albright's  course,  E.E.  355.   The  students  took  the 
lessons  on  PLATO  instead  of  attending  the  first  few  lectures  of  the  course. 
In  the  opinion  of  the  students,  the  machine  lessons  were  a  great  help  in  com- 
pleting their  learning  and  understanding  of  the  material.   The  material  in 
the  lessons,  however,  was  not  new  to  them.   The  real  test  of  learning,  i.e., 
subjecting  the  students  to  new  material  and  testing  their  comprehension,  could 
not  be  evaluated  in  this  experiment . 

E  .  R  .  Lyman 

5.5  PLATO  Laboratory 

A  computer  program  has  been  written  which  allows  the  students  to  make 
simulated  laboratory  measurements.   These  measurements  will  provide  necessary 
data  for  problems  posed  in  the  main  PLATO  text.   When  a  student  presses  a 
LAB  button,  he  goes  to  a  LAB  routine  just  as  he  ordinarily  would  go  to  HELP 
when  pressing  a  HELP  button,   In  the  LAB  he  can  construct  simple  experiments 
to  obtain  the  measurements  he  needs  .   The  computer  performs  the  necessary 
calculations  to  make  the  indicated  measurements  and  provides  the  student  with 
the  desired  data.   The  student  can  make  as  many  measurements  as  he  wishes 
until  he  has  sufficient  data  to  answer  the  problem.   He  then  returns  to  the 
main  PLATO  sequence  . 

The  program  has  been  written  as  a  general  subroutine  for  PLATO  II  so 
that  it  may  be  used  for  a  variety  of  simple  experiments  by  changing  only  the 
parameters  from  experiment  to  experiment . 

D.  L.  Bitzer 

E .  R .  Lyman 

5.6  Recording  of  Physiological  Measurements 

An  exploratory  study  was  launched  in  September,  1962,  concerning  physio- 
logical reactions  of  students  to  various  types  of  discovery-producing  problem 
sequences.   The  strategy  of  the  study  is  to  record  physiological  variables 
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such  as  heart  rate,  galvanic  skin  response,  and  parameters  of  the  breathing 
cycle  and  look  for  changes  which  occur  in  conjunction  with  discovery  exper- 
iences .   The  increased  motivation  which  seems  to  accompany  the  use  of  well- 
prepared  discovery  exercises  has  been  attributed  to  "intrinsic  reinforcement" 
occurring  in  concurrence  with  "discovery."   If  "intrinsic  reinforcement" 
events  do  occur,  it  should  be  possible  to  detect  and  record  their  physio- 
logical correlates .   Heart  rate  was  selected  as  the  first  physiological 
variable  for  study.   It  is  being  measured  as  the  reciprocal  of  the  interval 
between  two  consecutive  spikes  on  the  student's  ECG  trace. 

Instrumentation  for  simultaneous  recording  by  stop-motion  photography 
of  the  student's  heart  rate  and  of  the  textbook  and  blackboard  has  been 
developed.   Films  have  been  made  on  several  students  while  they  were  working 
on  ILLIAC  words.   A  "chapter"  on  number  sequences  was  prepared.   It  has  been 
tried  out  with  two  students  and  revised.   The  revised  "chapter"  and  tapes 
have  been  checked  on  one  student  and  corrected.   Thus,  we  are  now  ready  for 
data  collection  on  this  chapter.   Subsequent  chapters  will  be  designed  in 
order  to  investigate  the  different  emotional  response  of  students  to  different 
types  of  problem-sequences.   Analyses  of  preliminary  films  of  heart  rate  are 
being  conducted  to  determine  what  interesting  correlations  can  be  found  with 
events  in  the  PLATO  lesson. 

J.  A.  Easley,  Jr. 
R.  A.  Avner 

5.7  PLATO  III  Computer  Programs 

Work  continues  on  the  preparation  of  a  PLATO  III  program  for  the  CDC 
1604.   The  important  features  of  this  program  were  outlined  in  the  Quarterly 
Progress  Report  for  June  1962.   Currently,  programming  details  such  as  the 
organization  of  the  data  banks  are  being  worked  out.   At  the  same  time,  pro- 
grams are  being  prepared  corresponding  to  the  various  possible  student  inputs. 
It  is  expected  that  the  complete  PLATO  III  program  will  become  available  in 
the  next  quarter . 

P.  Braunfeld 
L.  Morgan 
M.  Secrest 

5.8  PLATO  III  System 

The  PLATO  III  system  is  a  computer  controlled  teaching  system  designed 
to  accomodate  up  to  32  student  sites  .   The  student  sites  are  connected  to 
one  input  and  one  output  channel  of  a  CDC  1604  digital  computer  by  means  of 
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an  electronic  interface.   The  PLATO  III  system  is  presently  in  the  process 
of  being  constructed,  with  plans  to  start  operation  with  two  students  early 
in  1963.   The  remainder  of  the  student  sites  will  be  added  as  they  are  com- 
pleted . 

5.8.1  Electronic  Interface 

Circuit  design  work  has  been  finished  and  circuit  panels  wired.   The 
panels  have  not  yet  been  assembled  as  modules  for  interconnection  and  testing. 
The  mechanical  design  and  construction  of  the  modules  is  nearing  completion 
and  will  soon  allow  assembly  of  the  modules. 

5.8.2  Keyset  Input  Logic 

Panels  have  been  wired,  ready  for  assembly  as  modules  and  testing. 

5.8.3  Output  Logic 

The  circuit  design  for  logical  control  of  the  storage  tubes  has  been 
finished.   Circuit  boards  for  one  student  are  almost  completed. 

Slide  select  and  special  effect  register  panels  for  two  sites  have  been 
wired,  but  not  assembled  as  modules.   The  two  video  switches  which  will  be 
operated  by  the  slide  select  registers  have  already  been  tested. 

5.8.4  Slide  Selector 

Work  on  the  second  slide  selector  has  been  temporarily  discontinued  in 
order  that  prototype  circuitry  for  the  PLATO  III  System  may  be  built.   The 
delay  thus  caused  in  completion  of  the  slide  selector  was  felt  justified  in 
that  there  appeared  no  pressing  need  for  a  second  slide  selector  at  the  pre- 
sent time. 

5.8.5  Storage  Tube  Circuitry 

During  this  quarter,  development  and  construction  of  circuitry  to  provide 

for  various  modes  of  operation  of  32  storage  tubes  utilized  by  the  PLATO  III 

system  continued.   Layout  and  fabrication  problems  regarding  video  pre-amplif ier, 

digital-analog,  and  deflection  amplifier  circuitry  were  considered  and  solved 

to  the  extent  that  construction  of  prototype  circuitry  was  begun.   It  is  expected 

that  the  aforementioned  prototype  circuitry  will  be  available  for  evaluation 

during  the  forthcoming  quarter. 

D.  Diamond 
B.  Voth 
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6  .   VACUUM  INSTRUMENTATION 


D.  Alpert  Jose  L.  de  Segovia 

W.  C.  Prothe  Harry  Tomaschke 

Donald  Lee  Nicholas  Vassos 
Wil  Schuemann 


6.1   Differential  Gauge  Studies 

With  the  development  of  the  inverted  ionization  gauge  by  R.  T.  Bayard 

22 
and  D.  Alpert   ,the  production  of  low  pressures  was  so  stimulated  that 

-10  23  24 

pressures  of  the  order  of  10    Torr  are  commonplace  today.   '     Further 

advances  in  attaining  even  lower  pressures  have  been  hampered  by  the  lack  of 

a  simply  constructed  gauge  which  is  capable  of  reading  pressures  below  the 

usually  accepted  x-ray  limit  of  the  Bayard-Alpert  gauge.   This  paper  describes 

a  modification  of  the  inverted  gauge  in  which  the  effective  x-ray  limit  has 

been  substantially  reduced. 

In  this  modif Ication,  a  second  collector  has  been  added  to  the  gauge 

and  is  biased  at  a  potential  below  that  of  the  first  collector  (see  Figure 

,25 

10)   .   With  a  bias  of  sixty  volts,  the  ion  current  from  the  biased  collector 

was  found  to  be  twice  that  of  the  unbiased  collector.   On  the  other  hand,  the 
x-ray  currents  from  these  collectors  are  approximately  equal  due  to  their 
similar  construction  and  symmetric  placement  of  the  collectors  within  the 
gauge.   Thus,  if  the  collector  currents  are  subtracted,  the  x-ray  currents 
nearly  cancel  and  the  net  current  is  to  first  order  due  only  to  ions.   Since 
the  ion  current  to  the  biased  collector  is  twice  that  due  to  the  unbiased 
electrode,  the  difference  current  corresponds  to  that  of  a  gauge  with  a 
sensitivity  one-third  that  of  a  standard  gauge,  but  with  a  greatly  reduced 
x-ray  current . 


22 

R.  T.  Bayard,  and  D.  Alpert,  Rev.  Sci  .  Inst.  2JL,  571  (1950). 

D.  Alpert,   Production  and  Measurement  of  Ultrahigh  Vacuum,"  Springer- 

Verlag  (Berlin),  1958,  p.  609. 

24 

H.  A.  Steinherz  and  P.  Redhead,  "Ultrahigh  Vacuum,"  Scientific  American 

206,  78  (1962)  . 
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Redhead  has  suggested  a  gauge  involving  the  use  of  a  fourth  element  to 

modulate  the  ion  current  with  an  ac  voltage.   P.  Redhead,  Rev.  Sci.  Inst. 

31,  943  (1960) . 
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A  modified  gauge  was  constructed  by  a  simple  modification  of  a  commercial 
version  of  the  Bayard-Alpert  gauge.   The  bias  voltage  typically  used  is  60 
volts,  though  this  is  not  critical.   The  design  of  the  biasing  circuit  is, 
however,  quite  important,  since  leakage  currents  can  be  as  large  as  the  ion 
currents  to  be  measured.   However,  it  has  been  found  that  with  the  proper 
use  of  guard  ring  potentials  on  the  gauge,  as  well  as  on  the  biasing  device, 
the  leakage  currents  may  be  reduced  to  tolerably  low  values.   Additionally, 
the  use  of  low  noise  coaxial  cable  was  found  to  be  necessary. 

A  Keithly  differential  electrometer  was  found  to  be  ideally  suited  to 
perform  the  subtraction  of  the  two  collector  currents. 

The  modified  gauge  was  tested  on  an  ultrahigh  vacuum  system.   At  a 

-11 
gauge  reading  corresponding  to  a  few  times  10    Torr,  the  net  current  was 

obtained  as  a  function  of  the  energy  of  the  ionizing  electrons.   Figure  11 
indicates  two  typical  runs.   The  curve  labeled  one  (1)  was  secured  by  biasing 
collector  one,  whereas  the  curve  labeled  two  (2)  was  secured  by  biasing  the 
second  collector  through  the  interchange  of  the  electrometer  leads  at  the 
ion  gauge.   Some  of  the  remaining  contribution  of  the  x-ray  current  is  evi- 
denced by  the  slowly  increasing  tail  in  curve  one  and  the  slowly  decreasing 
tail  in  curve  two.   The  interchange  of  the  leads  explains  the  change  in  sign 
of  the  contribution  of  the  x-ray  current.   Since  the  contributions  are  of 
opposite  signs,  however,  the  sum  of  curves  one  and  two  should  further  reduce 
the  x-ray  current.   Curve  three  (3)  is  this  sum  and  is  the  general  shape  of 
an  ionization  curve  up  to  about  400  volts  . 

The  addition  of  the  second  collector  allows  this  gauge  to  be  used  in 
the  Redhead  mode  of  operation;  that  is,  the  second  collector  may  be  used  as 
a  modulating  element.   In  this  mode  of  operation  the  x-ray  limit  is  approxi- 
mately the  same  as  in  the  differential  mode.   However,  the  differential  mode 
of  operation  is  free  of  troublesome  transients  introduced  from  the  necessary 
switching  or  modulation.   Furthermore,  the  differential  mode  provides  continu- 
ous pressure  readings  with  time  constants  of  a  few  seconds,  the  time  constant 
of  the  electrometer,  rather  than  some  multiple  of  the  switching  frequency 
which  is  roughly  2  cycles  per  minute. 
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Figure  11.   Collector  currents  in  double-collector  gauge 
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The  modified  gauge  has  a  relative  x-ray  limit  at  least  one,  and  possibly 
two,  orders  of  magnitude  below  a  standard  inverted  ionization  gauge.   The 
usefulness  of  the  modified  gauge  is  limited  by  the  quality  of  the  differen- 
tial micro-ammeter  and  by  the  ability  of  the  experimenter  to  eliminate  stray 
dc  leakage  currents .   Since  the  gauge  is  a  simple  modification  of  the  original 
Bayard-Alpert  design,  all  power  supplies  remain  suitable.   Further,  since 
quality  differential  micro-ammeters  are  available  commercially,  this  guage 

represents  an  attractive  solution  to  the  problem  of  reading  pressures  into 

-12 
the  10    Torr  pressure  region. 

The  author  wishes  to  acknowledge  the  many  valuable  discussions  with  Dr. 
Daniel  Alpert,  whose  comments  immeasurably  helped  in  the  development  of  this 
gauge  . 

D.  Lee 
6.2   Vacuum  Gauge  Studies 

Progress  in  the  development  of  the  photo  current  suppressor  gauge  for 
the  measurement  of  total  pressures  in  vacuum  systems  has  been  reported  in 
three  previous  progress  reports  .  The  development  of  the  gauge  has  for  the 
present  been  finished.  We  are  in  the  process  of  building  ten  of  the  model 
18  gauges  for  use  at  C.S.L.  Five  gauges  have  been  completed  and  are  being 
used  on  vacuum  systems . 

Since  the  presentation  of  a  paper  on  this  gauge  at  the  American  Vacuum 
Society  conference,  fifty  reports  along  with  detailed  construction  drawings 
have  been  sent  to  workers  in  the  vacuum  field  who  have  requested  them.   It 
seems  that  the  photo  current  suppressor  gauge  shows  much  promise  of  being  a 
useful  scientific  research  tool.   The  gauge  is  shown  in  Figure  12. 

Recently  several  papers  have  been  published  which  report  anomalous  be- 
havior in  the  Bayard-Alpert  inverted  ionization  gauges  when  operated  at  low 
pressures,  low  emission  currents,  or  when  exposed  to  certain  gases.   We  are 
trying  to  duplicate  the  published  results  and  plan  to  look  for  the  explana- 
tion of  the  phenomena  responsible.   Toward  this  end,  one  system  has  been 
modified,  and  a  second  system  has  been  built  and  is  now  operating.   Prelimi- 
nary data  shows  that  there  is  indeed  some  anomalous  behavior  concerning  the 
sensitivity  and  the  x-ray  limit  of  Bayard-Alpert  gauges  at  residual  pressures 
in  the  10    Torr  region. 

Wil  Schuemann 


47 


Figure  12  .   Model  18  photo  current 


suppressor  gauge 
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6.3   Total  and  Partial  Pressures:   Mass  Spectrometer  Studies 

Comparative  studies  of  total  and  partial  pressures,  using  two  types  of 
mass  spectrometers  and  a  Bayard-Alpert  gauge,  have  continued  during  this 
quarter.   The  ultrahigh  vacuum  system  used  is  not  entirely  suitable  for  ac- 
curate quantitative  results  because  of  several  small  conductances.   However, 
some  interesting  data  have  been  obtained. 

6.3.1   The  Ion  Currents  for  Ar,  Ne,  and  N0,  as  Indicated  by  Both  Mass 

Spectrometers  and  Bayard-Alpert  Gauge,  Have  Been  Studied 

The  following  points  can  be  made: 

(a)  Relations  of  the  form 

+        +   "l      + 
I  ,  -   a, (I„.)   ,     I„„  =  ion  current  as  indicated  by  Davis-Vanderslice, 

DV     1   (3A    '      DV  ' 

+       +2       + 
In  =  a  (I   )   ,      In  =  ion  current  as  indicated  by  Omegatron, 

o  l     2   pA  db 

I   =  ion  current  as  indicated  by  Bayard-Alpert  gauge 

PA 

have  been  found  in  which  the  exponents  |3   and  (3   are  not  unity.   These  expo- 
nents vary  with  the  electron  current  of  the  Bayard-Alpert  gauge  and  with 
whether  one  of  the  mass  spectrometers  was  turned  on  or  off. 

The  exponent  value  of  unity  lies  outside  of  the  probable  experimental 
error,  except  in  one  case.   This  case  consisted  of  introducing  an  Argon  leak 
into  the  system  with  the  Davis-Vanderslice  turned  off  .   For  those  gases  studied 
the  exponents  ranged  between  0.9  and  1.1. 

(b)  In  the  relations  of  the  form 

the  exponent  (3   is  closer  to  unity  than  in  the  last  case. 

(c)  When  pure  N  was  leaked  into  the  system,  an  appreciable  increase 
of  peaks  16  and  12  was  noted  by  both  spectrometers.   This  means  that  the 
partial  pressure  of  CO  in  the  system  was  increased. 
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(d)   It  was  also  observed,  if  we  assume  that  the  ion  current  is  pro- 
portional to  the  pressure  in  each  spectrometer,  that  the  cracking  patterns 
studied  are  not  linear  with  the  pressure. 

6.3.2  Residual  Vacuum 

(a)  The  residual  vacuum  indicated  by  each  spectrometer  was  a  function 
of  the  spectrometer  itself. 

(b)  The  influence  of  both  spectrometers  on  the  residual  vacuum  has  been 
studied.   The  Davis-Vanderslice  has  been  found  to  have  a  more  important  con- 
tribution to  the  peak  28  than  the  Omegatron. 

6.3.3  Recent  Experimenting 

It  was  also  mentioned  in  the  previous  Progress  Report  that  the  Davis- 
Vanderslice  mass  spectrometer  revealed  an  anomalous  mass  16.   During  these 
experiments,  more  information  about  this  anomalous  peak  was  sought.   The 
following  considerations  were  revealed  by  experiment: 

(a)  An  equilibrium  between  the  partial  pressure  of  CO  and  the  ion 
current  for  mass  16  exists. 

(b)  Outgassing  the  grid  of  the  Davis-Vanderslice  with  electrons  of  high 


energy  (1200  eV,  1.5  mA)  did  not  reveal  any  new  information. 

167  28 
on  the  conditions  of  the  system  and  the  spectrometer. 


(c)   The  relation  Inc/I00  has  been  found  to  be  from  0.5  to  0.15,  depending 


(d)  When  the  partial  pressure  of  CO  in  the  system  was  increased,  the 

+     4-  +     + 

relations  I,_/I00  and  I,0/Ioo  were  found  to  be  reasonable. 
16   28       12   28 

(e)  If  the  pressure  of  air  introduced  into  the  system  is  maintained  at 
1  X  10   Torr,  during  the  first  minutes  all  peaks  12,  14,  16,  28  increased; 
after  this  period,  peak  12  decreased  and  16  continued  increasing.   The  time 
constant  for  this  increase  was  very  long  (on  the  order  of  two  hours)  . 

(f)  It  seems  reasonable  to  assume  that  this  phenomenon  occurs  on  the 
grid  of  the  Davis-Vanderslice  mass  spectrometer.   However,  more  data  and  a 
better  knowledge  of  the  conditions  of  the  system  are  necessary  to  explain 
the  phanomenon . 

All  these  results  are  being  consolidated  and  some  experiments  will  be 
repeated . 

Jose  L.  de  Segovia 


Also  at  Instituto  'L,  Torres  Quevedo,"  Madrid,  Spain. 
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C  .  Chen 

J  .  Crowder 

J.  Harris 

B.  Hicks 


T.  Lie 

C.  Mendel 

H.  G.  Slottow 

J.  B.  Tsui 


7.1   Elastic  Collisions  of  Electrons  with  Excited  Atoms 

In  the  course  of  assembling  the  instrumentation  for  the  experiment,  a 

Polarity  Coincidence  Detector  was  constructed  and  tested.   The  theory  of 

26 

operation  has  been  described  partly  in  the  literature.     Additional  theo- 
retical considerations  will  be  presented  in  a  forthcoming  report. 

The  schematic  diagram  is  sketched  in  Figure  13.   The  signal  of  frequency 

oj  plus  the  superimposed  noise  is  passed  through  a  limiting  amplifier  which 
o 

converts  it  into  a  square  wave.   The  reference  signal  is  converted  into  narrow 
pulses  whose  repetition  frequency  is  equal  to  the  signal  frequency.   One  in- 
phase  and  one  quadrature-phase  channel  are  provided.   The  limited  signal-plus- 
noise  is  sampled  with  these  pulses,  and  coincidences  are  recorded  whenever 
the  polarity  of  the  wave  form  under  investigation  coincides  with  the  polarity 
of  the  sampling  pulse. 

With  narrow  band  gaussian  noise  present,  the  fraction  of  sampling  pulses 
resulting  in  a  coincidence  is  given  by 


n(0)  =  i  [1  -  erf  (x   -  A// 2  cos  <j>)/0"  /[/  2] 


for  the  in-phase  channel,  and 


n(7T/2)  =  |[  1  -  erf  (x^  -  A  //  2  sin  $)/<J   //  2  ] 


for  the  quadrature-phase  channel.   The  symbols  are  A  for  the  r.m.s.  signal 
amplitude,  0"  for  the  variance  of  the  noise,  and  x   for  the  discrimination 
level  (the  level  which  the  signal  has  to  exceed  in  order  to  be  registered 
as  a  coincidence) . 
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Ben  S.  Melton  and  P.  R.  Karr,  Geophysics  22,  (1957),  p.  553 
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For  small  values  of  the  argument  (small  signal-to-noise  ratio  and  hard 
limiting)  the  error  function  can  be  expanded  and  the  response  of  the  detector 
becomes  linear.   The  signal-to-noise  ratio  is  then  given  by 


f 


2  2\   1/2 

[n(0)  -  nN(0)]   +  [n(77/2)  -  1^(77/2)]  f      =  A/(a/  IT), 


in  which  n   is  the  fractional  count  with  only  noise  present.   Figure  14  shows 
a  plot  of  the  measured  output  versus  input  signal-to-noise  ratio  (S/N) .   The 
signal  frequency  was  10  Kc/sec. 

Of  interest  is  the  statistical  error  of  the  measurement  .   If  the  samples 
were  uncorrelated,  the  statistical  error  of  N  counts  would  be  AN  =  v  N.   If 
the  noise  is  passed  through  a  narrow  band  filter  prior  to  detection,  however, 
the  samples  are  correlated  and  the  statistical  error  is  increased.  We  write 

AN  =  p  VN, 

where  (3  is  a  function  of  the  band  width  and  the  form  of  the  spectral  density. 
The  experimental  values  are  compared  with  two  models 

-,2 
1)  S(gj)  =  S   exp  -  [  (u  -  u  )/Aw]   , 

o  o 


2)      S(cj)  =  S  \  1  +  c[  (go  -  oo  )/Aw]  (   , 
°  I  °       J 


c  =  e  -  1  =  1.71 


The  factor  e  -  l  has  been  added  to  normalize  the  significance  of  Aw  as  the 
1/e  -  width. 

In  Figure  15  measured  experimental  values  are  compared  with  these  two 
models.   The  deviation  is  apparently  due  to  the  insufficient  knowledge  of  the 
actual  spectral  density.   Qualitatively,  however,  it  is  demonstrated  that  the 
statistical  error  increases  with  decreasing  band  widths. 

Summarizing,  we  can  say  that  the  developed  Polarity  Coincidence  Detector 
works  reliably  for  signals  in  the  range  20  c/sec  to  20  Kc/sec  and  performs  ac- 
cording to  theoretical  expectations. 

M.  Raether 
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Figure  14.   Detection  curve  for  polarity  coincidence 
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Figure  15.   Comparison  of  experimental  data  with 

models  1  and  2  of  polarity  coincidence 
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7.2   Atomic  Processes  in  Noble  Gas  Mixtures 

In  this  quarter,  measurement  of  electron-xenon  collision  probability 
for  momentum  transfer,  P  (Xe),  at  low  energies  is  completed.   The  method  and 
the  subsequent  treatment  of  the  result  of  the  experiment  have  been  reported 
in  the  previous  report  and  will  not  be  repeated  here. 

The  measured  values  of  the  effective  electron  collision  frequency  per 

unit  xenon  gas  pressure,  v(Xe)/P  (Xe),  as  a  function  of  electron  temperature 

o 

is  shown  in  Figure  16.   Adopting  the  same  functional  form  of  the  cross  section 
on  electron  velocity  as  in  the  case  for  krypton,  the  best  fit  to  our  experi- 
mental data  yields 


v(Xe\    =    5.59  x  108T  V2  -  3.69  x  10  V  +  6.96  X  105T  3/2       (7.1) 
p  (Xe)  e  3  e 


and  is  shown  as  the  solid  curve  on  Figure  16  .   The  energy  dependence  of  P 

m 

calculated  from  Equation  (7.1)  and  plotted  conventionally  is  shown  in  Figure  17, 
Data  obtained  by  C.  Ramsauer  and  by  C.  Ramsauer  and  R.  Kollath  at  higher 
energies  are  also  shown. 

A  detailed  report  on  atomic  processes  in  HeKr  and  HeXe  mixtures  is  under 
preparation.   The  same  technique  is  applied  to  other  cases  in  studying  low 
energy  electron  and  ion  behavior  in  plasmas . 


Two  new  projects  have  been  initiated  in  this  quarter;  namely,  the  forma- 
process  of  N    through  two  body  collisj 
propagation  of  ionic  sound  waves  in  plasmas 


+  +  + 

tion  process  of  N    through  two  body  collisions  N    +  N   ^N,   and  the 

*  4  2     2        4 


C.  L.  Chen 
7.3   The  Boltzmann  Equation  for  a  Shock  Wave 

The  first  complete  1604  program  has  been  in  use  during  the  past  quarter. 
It  will  be  recalled  that,  as  compared  to  the  ILLIAC  program,  this  program  had 
been  improved  in  two  ways:   by  more  accurate  moment  calculation  and  by  correc- 
tion of  an  incorrect  mask.   Many  calculations  have  now  been  made  for  a  Mach 
number  of  3.17,  for  J  =  9,  and  17  stations  in  the  shock  wave.   For  J  =  9, 
the  principal  results  are  essentially  the  same  as  for  the  ILLIAC  runs.   For 
J  =  17,  the  iterations  now  appear  to  be  stable,  though  they  were  unstable  for 
the  corresponding  ILLIAC  runs  .   Thus  apparently  one  of  the  questionable  features 
of  the  ILLIAC  results  has  been  eliminated. 
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The  1604  calculations  are  faster  by  about  a  factor  of  15  as  compared  to 
corresponding  ILLIAC  calcualtions.   This  greater  speed  makes  possible  many 
types  of  study  that  were  not  feasible  before . 

One  of  several  new  and  useful  auxiliary  programs  will  be  mentioned. 
This  program  is  written  in  FORTRAN  language,  and  it  forms  means  and  standard 
deviations  of  the  moments  of  the  velocity  distribution  functions  in  the  out- 
put of  the  1604  machine  language  program.   It  also  calculates  the   t-test  . 
Statistical  analysis  of  the  iterations  can  now  be  made  in  a  few  minutes  that 
required  several  months  of  hand  calculations . 

In  the  Department  of  Aeronautical  and  Astronautical  Engineering,  Professor 
S.  M.  Yen  and  a  graduate  student  are  starting  a  detailed  comparison  of  the  many 
different  calculations  of  shock  structures  that  have  been  made  by  various 
investigators  in  the  past  15  years.   Their  studies,  though  not  supported  directly 
by  CSL,  are  being  coordinated  with  the  CSL  1604  calculations,  with  mutual  bene- 
fit to  both  groups  , 

In  the  next  quarter,  we  expect  to  use  the  1604  program  for  calculations 
for  lower  Mach  numbers,  in  order  to  find  if  the  results  here  are  reasonable 
as  they  were  not  for  the  ILLIAC  program,  and  to  provide  a  basis  for  direct 
comparison  with  experiment  and  with  analytical  theories  of  shock  structure. 
We  also  plan  to  start  the  study  of  the  accuracy  of  the  Monte  Carlo  calcula- 
tions of  the  collision  integral  as  a  function  of  sample  size,  this  study  being 
feasible  with  the  much  greater  speed  of  the  1604  program. 

B.  Hicks 

J.  K.  Aggarwal 

7.4   Plasma  Delay  Line 

The  drift  of  charged  particles  in  crossed  electric  and  magnetic  fields 
can  in  principle  provide  electronically  controllable  delays  of  information 
for  wide  band  communication  systems.   Two  approaches  to  design  are  being 
studied.   In  one,  electrons  and  ions  neutralize  the  space  charge  as  they 
drift  to  the  collector  region.   In  a  second,  electrons  alone  carry  the 
delayed  signals.   Some  results  of  this  investigation  have  been  discussed  in 
previous  CSL  Progress  Reports .   In  this  section  we  consider  electric  field 
structures  appropriate  for  an  electron  drift  line.   We  assume  that  the  elec- 
tron beam  is  continuous,  that  the  drift  space  contains  no  conductors,  and  that, 
to  avoid  gradient  drifts,  the  magnetic  field  is  constant  in  the  drift  space. 


59 


The  first  requirement  for  an  electron  drift  delay  line  is  that,  for  all 
electrons,  the  times  of  flight  from  injector  to  collector  be  substantially 
equal.   Although  in  rectangular  geometry  this  is  equivalent  to  the  require- 
ment that  the  drift  speeds  be  equal,  in  circular  geometry  these  requirements 
are  not  the  same,  and  it  is  the  implications  of  the  differences  that  interest 
us.   For  rectangular  drift,  constant  speed  implies  a  constant  transverse  com- 
ponent of  electric  field,  E  ,  which  in  turn  implies  that  E  ,  the  longitudinal 
component  of  electric  field,  i.e.,  parallel  to  the  magnetic  field,  must  be 
zero.   Since  the  electrons  are  not  focussed  in  this  direction  they  escape  from 
the  drift  space  along  the  magnetic  flux  lines.   For  circular  drift  the  fields 
can  be  arranged  to  provide  both  equal  times  of  flight,  and  also  focussing 
forces  along  the  magnetic  field  lines. 

Laplace's  equation  is  cylindrical  coordinates,  when  the  angular  depen- 
dance  is  zero,  has  the  form: 

3e     E     9E 
r     r      z 

-9T  +  —  +  ~3z   =   °> 


r 


where 

E  -   radial  component  of  the  electric  field, 
E   =  longitudinal  component  of  the  electric  field, 
r  =  radial  coordinate, 
z  =  longitudinal  coordinate. 
The  magnetic  field  B  is  assumed  to  be  directed  along  the  z  axis,  and  the 
plane  of  symmetry  toward  which  electrons  are  focussed  by  E   is  defined  by 
z  =  0.   It  is  easy  to  verify  that  a  solution  is 


E   =  -ar 
r 

E   =  2az 
z 

If  a  >  0,  E   points  away  from  the  plane  z  =  0,  and  electrons  will  be  focussed. 
The  drift  velocity  of  the  electrons,  v  =  E  /B  =  E  /B,  increases  linearly  with 
the  radius.   Since  the  electrons  are  tied  to  equipotential  lines  by  the  magnetic 
field,  this  is  precisely  the  variation  needed  to  maintain  equal  times  of 
flight  from  an  injector  at  one  angle  to  a  collector  at  another. 
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A  corresponding  equation  for  the  potential  is 


V(r,z)  -  V(r  ,0)  =  (a/2)(r2  -  r  2)  -  az2, 
o  o        ' 


where 

V(r,z)  =  potential  at  the  point  (r,z),  and 

r   =  an  arbitrary  radius . 

o  J 

The  equipotential  surfaces  are  therefore  hyperboloids  of  revolution.   A 
suitable  structure  for  a  delay  line  could  be  formed,  Figure  18,  from  sections 
of  two  hyperboloids  together  with  field  terminations  in  the  planes  of  the 
section  edges  . 

We  have  so  far  neglected  the  effects  of  space  charge.   If  we  now  assume 
that  the  electron  beam  is  solid  with  a  uniform  charge  density,  we  can  obtain 
simple  expressions  for  the  electric  fields  and  the  potentials.   Poisson's 
equation  is,  for  this  case, 

9e     E     3e 

r     r      z   P    ,, 


where 

P  =  change  density 


6  =  dielectric  constant, 

-12 
=  8.85  x  10    farads/meter, 


The  radial  electric  field  again  has  the  form 


E   =  -ar, 
r 


but  the  longitudinal  field  has  an  additional  term: 


E   =  (2a  +  P/e)z 
z 
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The  equation  for  the  potential  is  now 


V(r,z)  -  V(r  ,0)  =  (a/2)(r2  -  r  2)  -  (a  +  p/2e)z2 
o  o 


The  equipotential  surfaces  are  again  hyperboloids  of  revolution,  but  as  the 
magnitude  of  the  charge  density  increases  (p  more  negative),  the  angle  be- 
tween the  asymptotic  cones  increases  until  -p/2e  =  a.   For  this  density, 

the  equipotentials  are  the  circular  cylinders,  E   =  0,  and  the  longitudinal 

z 

focusing  disappears . 

To  obtain  some  idea  of  the  magnitudes  of  the  variables  in  a  realistic 
case,  we  consider  the  following  parameters,  with  -  p/E   =  a: 

r  =  0.1  meter, 

-12  3 

P  =  10    coulombs/cm  , 

B  =  5, 000  gauss . 

We  find  that 

5 
E   =  0.113  x  10  volts/meter, 
r  7 

-6  2 

(Current  density)  =  2.26  x  10   amperes/cm  , 

(time  delay  for  one  circular  excursion)  =  27.8  x  10   seconds. 
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8.   HIGH  MAGNETIC  FIELD  SUPERCONDUCTORS 


C.  B.  Satterthwaite  W.  C.  Leung 

J.  0.  Kopplin  R.  Ries 

E.  Wisman 


8.1  Introduction 

This  group  is  primarily  interested  in  studies  of  superconductivity 
pertaining  to  the  persistance  of  superconductivity  in  very  high  magnetic 
fields  .   There  are  three  phases  of  this  program  in  which  work  is  being 
initiated.   They  are:   (a)  basic  studies  of  superconducting  phenomena, 
(b)  preparation  and  study  of  materials,  and  (c)  devices  and  apparatus  involving 
superconductors  . 

8.2  Superconducting  Solenoid 

The  superconducting  solenoid  described  in  the  preceding  progress  report 

was  subjected  to  further  testing.   It  was  found  that  resistance  in  the  joints 

between  sections  of  the  windings  prevented  reaching  the  expected  magnetic 

field  for  a  given  current  .   With  resistance  in  the  superconducting  path  part 

of  the  current  flowed  through  the  copper.   Some  improvement  has  been  made 

by  making  superconducting  joints;  fields  in  excess  of  25,000  gausses  have 

been  generated.   It  appears,  however,  that  a  solenoid  of  this  design  will 

not  be  a  very  useful  laboratory  tool  because  of  the  very  large  time  constant . 

C.  B.  Satterthwaite 
J.  0.  Kopplin 
W.  C.  Leung 
R.  Ries 

8.3  Preparation  of  Nb  Sn 

Some  preliminary  experiments  indicated  that  one  can  get  crystalline 
compounds  of  Nb  and  Sn  by  controlled  cooling  of  a  solution  of  Nb  in  molten  . 
Sn  at  high  temperatures  or  by  maintaining  it  in  a  temperature  gradient.   A 
furnace  has  been  built  somewhat  similar  to  a  silicon  crystal  furnace  for 
maintaining  a  saturated  solution  of  Nb  in  Sn  at  1100  to  1200  centigrade 
degrees,  into  which  can  be  lowered  an  Nb  probe  which  will  be  at  a  somewhat 
lower  temperature  and  presumably  will  initiate  crystal  growth.   It  is  hoped 
that  single  crystals  of  Nb  Sn  and  perhaps  other  Nb-Sn  compounds  can  be  pro- 
duced in  this  way  for  studies  of  their  superconducting  properties. 

W .  C .  Leung 
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8.4  Vanadium  Study 

A  program  has  been  initiated  to  study  the  superconducting  properties 
of  vanadium.   Vanadium  was  selected  because  it  is  a  member  of  the  group  of 
"hard"  superconductors,  can  presently  be  obtained  in  relatively  pure  form 
and  holds  promise  of  having  interesting  superconducting  properties . 

A  quarter  inch  diameter  vanadium  rod  was  passed  through  an  electron 
beam  zone  refining  operation  several  times  for  impurity  removal.   The  speci- 
men was  then  swaged  back  into  a  true  cylindrical  shape,  and  a  single  crystal 
approximately  one  fifth  or  an  inch  in  diameter  and  3  1/2  inches  long  was 
grown  using  the  electron  beam  zone  melting  method.   This  crystal  was  cut  into 
two  unequal  sections  by  the  spark-cutting  technique.   Back  reflection  x-ray 
patterns  made  from  the  cut  surface  indicate  a  single  crystal  and  a  negligible 
amount  of  surface  damage  due  to  cutting.   Present  plans  call  for  drawing  the 
smaller  section  to  determine  the  rate  of  hardening  due  to  cold  work  and  for 
cutting  two  parallel  wires,  approximately  two  inches  in  length,  from  the 
larger  section  to  be  used  for  future  measurements  of  the  superconducting 
properties  . 

J.  0.  Kopplin 

8.5  Niobium  Study 

Work  has  been  started  to  develop  a  reliable  procedure  for  preparing  the 
niobium  specimens  suitable  for  electron  transmission  microscopy  studies. 
This  program  has  been  undertaken  to  determine  the  feasibility  of  obtaining 
information  from  electron  transmission  microscopy  studies  which  would  be 
useful  for  relating  the  changes  in  the  superconducting  properties  with  the 
changes  in  the  material  structure  brought  about  by  cold  working  of  the  ma- 
terial.  The  initial  efforts  has  been  directed  toward  developing  satisfactory 
etching  and  masking  techniques  to  be  used  in  preparing  specimens  of  1000  A 
or  less  in  thickness  . 

J  .  0  .  Kopplin 

8.6  Liquid  Helium  Level  Indicator 

A  patent  application  has  been  prepared  for  the  highly  satisfactory 

liquid  helium  indicator  described  in  the  last  progress  report . 

R.  Riess 

C.  B.  Satterthwaite 
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8.7   Critical  Current  Measurements 

A  superconducting  materials  program  in  cooperation  with  the  Fansteel 

Metallurgical  Corp.  has  been  in  progress  for  some  time.  In  this  program 

Fansteel  has  agreed  to  undertake  preparation  of  certain  high  field  super- 
conductors and  the  cryogenic  group  in  CSL  will  measure  critical  currents 

and  critical  fields.   Some  materials  have  been  produced  and  a  program  of 
making  critical  current  studies  is  in  progress. 

C.  B.  Satterthwaite 
R.  Ries 
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VACUUM  ELECTRICAL  BREAKDOWN 


E .  M .  Lyman  F .  K .  Konrad 

T.  C.  Casale  H.  E.  Tomaske 


Electron  microscope  studies  of  electrode  surfaces  have  been  continued 
with  a  microscope  of  recently  improved  resolution.  *    .   In  order  to  obtain 
reproducible  data  it  was  necessary  to  perform  the  electrical  measurements  in 
a  "clean  system;"  i.e.,  a  well  out-gassed  cathode  in  an  ultrahigh  vacuum. 

To  do  this,  a  cathode  assembly  was  constructed  which  can  be  mounted  on 
a  flange  and  inserted  into  the  ultrahigh  vacuum  system.   After  the  electrical 
measurements  have  been  completed,  the  cathode  and  assembly  can  be  removed 
and  placed  in  a  modified  specimen  holder  for  viewing  in  the  electron  micro- 
scope.  The  cathode  itself  is  the  tip  of  a  .003  inch  tungsten  wire,  the  other 
end  of  which  is  spotwelded  to  a  filament.   The  anode  is  the  electropolished 
end  of  an  0.1  inch  tungsten  rod. 

Some  of  the  results  obtained  are  as  follows: 

(1)  A  weld  between  the  anode  and  cathode  due  to  physical  contact  in 

vacuum  (for  gap  measurement  purposes)  results  in  sharp  points  when  the  elec- 

-5 
trodes  are  separated.   Some  of  these  points  were  about  2-5  x  10   cm.  high 

-5 
with  a  tip  diameter  of  10    cm. 

(2)  Heating  the  cathode  to  about  2000  C.  smooths  the  tip  considerably, 

-5 
leaving  more  nearly  spherical  projections  with  diameters  of  6  x  10    cm.  and 

larger  . 

(3)  Prebreakdown  pulses  were  observed  with  pulse  lengths  varying  from 

-1      -5  ,  -3 

10    to  10   sec.   For  a  typical  case  (gap  spacing  d  =  7.5  x  10    cm.;  applied 

voltage  V  =  4  x  10  volts),  the  average  (dc)  current  was  5  x  10    amp.  with 

-6 
pulses  corresponding  to  2  x  10   amp.   Although  there  were  no  visible  discharges, 

operating  the  electrodes  under  these  conditions  for  15  minutes  resulted  in  a 


change  of  shape  of  the  cathode  corresponding  in  places  to  a  removal  of  a  layer 

-4 
of  material  about  10    cm.  in  thickness  or  approximately  3000  atomic  layers. 

(4)   Preliminary  measurements  of  current-voltage  relationships  indicate 

2 
that  plots  of  log  (I/V  )  vs  I/V  are  straight  lines.   Calculation  of  the 

field  enhancement  factors,  B,  (effective  field  E  „„   =  B  V/d)  from  the 

'  ef  f  . 


*RCA  EMU-2E 
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slopes  of  these  plots  gives  values  ranging  from  30  to  130,  depending  on 
the  previous  conditioning  of  the  cathode.   This  effect  will  be  investigated 
more  thoroughly. 

H.  E.  Tomaschke 
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10.   PHYSICAL  ELECTRONICS 


R.  N.  Peacock  G.  Riddle 


Evaporated  Thin  Films 

Thin  film  sandwich  structures  of  Al-Al  o  -Au  have  been  prepared  on 

glass  substrates  which  exhibit  negative  resistance  characteristics,  and  from 

which  vacuum  emission  of  hot  electrons  has  been  observed.   With  the  A1„0„ 

2  3 

layer  anodized  to  about  300  A  thickness,  electron  emission  currents  up  to 
one  microampere  can  be  obtained  from  a  2x2  mm  area.   The  voltage  between 
the  Al  and  Au  films,  for  emission,  is  very  near  the  breakdown  voltage,  and 
is  usually  in  the  range  of  10  to  15  volts,  for  an  oxide  film  a  few  hundred 
Angstroms  thick. 

The  diode  and  emission  currents  are  "noisy"  .   Fusing  of  the  thin  Au 
layer  has  also  been  a  problem. 

R.  N.  Peacock 


